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ABSTRACT 


A second-generation data acquisition and reduction syston is now being used 
at JPL to conduct long-baseline radio InterferMietry aeasurements. This system, 
called the Block 0 system, improves upon the previous 48-kbs system in many ways. 
Significant improvements Include a higher recording rate of 4 Hbits/s, multi- 
channel recording, multistation cross-correlation, phase calibration, S/X cali* 
bration of charged particle effects, and single-observation bandwidth synthesis. 
This report updates and improves the analysis previously performed for the 
48-kbs system and traces the form aiwl flow of information from signal reception 
to recorded bits at each station, through cross-correlation, phase-tracking, 
phase calibration, and bandwidth synthesis, to final observable formation. In 
this process, the data volume is greatly compressed for each observation, from 
'V'lO^ bits on the station tapes to only four final observables: amplitude, delay, 

phase-delay rate and RF ptv*se. It is shown that, after phase calibration and 
propagation media calibrations, the final delay observable is the sum of a geo- 
metric delay and a clock synchronization offset. 
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SECTION I 


INTRODUCTION 


Over the last few years, a second generation system for the acquisition and 
reduction of very long baseline Interferometry (VLSI) data has been developed 
at JPL. In most respects, this second generation digital system, the Block 0 
system. Improves considerably on the previous 48-kbs digital system, most signi- 
ficantly by providing (a) a higher record rate of 4 Mbit/s, (b) multichannel 
recording, (c) multistation cross-correlation, (d) phase calibration, (e) dual- 
band calibration of charged particle effects and (f) single-observation bandwidth 
synthesis. Because of these improvements, many changes have been required in the 
data reduction theory. This report augments and modifies previous analyses 
(Ret. 1) to encompass the new features. Even though many improvements have been 
Incorporated, some features of the BI.KO system were not optimally designed due 
to various circumstances that prevailed during the development period. A new 
system now under development, the Block II system, will greatly improve upon the 
BLKO system, more nearly realizing the potential of VLBI. 

The report Is organized as follows. Section II provides a summary of the 
data-processlng steps found in the BLKO system. Detailed derivation and analysis 
of those steps is presented in the subsequent sections. Section III gives 
fundamental definitions for basic time-keeping quantities. Section IV derives 
an expression for the recorded signal by tracing the signal from its free-space 
form through the instrumentation. In Section V, a non-relativistic model is 
developed for the "fast" fringes (cross-correlation function) produced by the BLKO 
correlator, the CIT/JPL Mark II correlator. In some respects, this part of the 
analysis repeats earlier work [1], but differs in important ways. Steps that are 
essentially the same are only summarized. One new feature is that instrumental 
terms are now grouped and defined in a more precise manner in order to demonstrate 
phase calibration through the use of tone phase. Further, the analysis treats 
delay offsetting of both bitstreams and provides an algorithm for correcting for 
the fractional-bit-shift errors introduced by the correlator. Section VI analyzes 
the fringe-stopping (counter-rotation) procedure used by the processor to com- 
press the data. Ii^,^^ection VII, the stopped fringes are converted by a Fourier 
transform from the iag^omain to the frequency domain. Transformation errors 
resulting from limited-lag reduction are discussed. Section VIII discusses phase 
calibration and derives an expression for the phase of a stopped tone. In Sec- 
tion IX, tone phase is used to correct fringe phase, thereby removing most instru- 
mental effects and casting the phase into a form subject to simple interpretation. 
In Section X, the fringe model used in phase tracking is formulated to match the 
frequency-domain fringes generated by the Fourier transform. As in the lag domain, 
there are four estimated parameters: fringe amplitude, bit-stream-alignment (BSA) 

delay, fringe phase and phase-delay rate. Output amplitude is absolutely normal- 
ized to first approximation by removing known constant amplitude factors. 

Further, it is shown that BSA delay and phase delay are cast in the same form 
with regard to instrumental effects. Finally, Sections XI-XIV present detailed 
descriptions of the steps that produce the final values for delay and phase-delay 
rate. Those steps include observable formation, model Improvement, propagation- 
media calibrations, and model restoration. It is shown that the final delay 
observable is the sum of a geometric delay and a clock synchronization offset. 
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SECTION II 

SUMMARY OF THE BLKO SYSTEM 


The BLKO system consists of several subsystems developed by a number of 
individuals from the Jet Propulsion Laboratory and the California Institute of 
Technology. The major subsystems and the people who developed them are data 
acquisition Instrumentation (E. J. Cohen, G. S. Parks, E. H. Slgman, L. J. 

Skjerve, and D. J. Spitzmesser) , phase calibrator (G. F. Lutes), correlator 
(M. S. Ewing and D. H. Rogstad) and postcorrelation software (G. H. Purcell, 

J. A. ScHeid and J. B. Thomas). This report presents a theoretical analysis of 
the ^Locessing of a random broadband signal from a natural source and traces 
the signal from its free-space form to the final output observables of amplitude, 
delny, and phase-delay rate. 

As shown schematically in Fig. 1, the random signal of a natural source 
received at an antenna at RF passes through various filters, amplifiers and 
mixers and a portion of it emerges at baseband spanning a 2-MHz-wlde passband. 
Details of the instrumentation are not discussed here since the Instrumental 
effects can be represented theoretically by a system bandpass function, a compo- 
site heterodyne signal, and composite phase shifts and delays. The instrumenta- 
tion can switch a heterodyne signal through a sequence of frequencies so that the 
baseband signal cycles through a number of channels corresponding to 2-MHz win- 
dows at selected frequencies at RF. In this "time-multiplexed" multichannel 
mode, the dwell time on each channel can be as short as 0.2 seconds and as long 
as several seconds. The baseband signal is infinitely clipped, sampled at a 
4-Mbit/s rate in a two-level mode and digitally recorded in a time-multiplexed 
format on magnetic tape by an IVC 825 or RCA 201 recorder. The several channels 
separated in frequency permit the subsequent reconstruction of a very wide 
effective bandwidth by a process called bandwidth synthesis (BWS). 

When tapes from the antennas are sent to the correlator for data reduction, 
they are played back at the original data rate and corresponding frequency chan- 
nels are cross-correlated. Multiple baselines (-10) and multiple channels (%20) 
can be simultaneously correlated. The output of the processor, called the stopped 
(or counter-rotated) fringes, is passed to post-correlation software that (a) 
transforms the fringes from the lag domain to the frequency domain, and (b) phase- 
tracks the transformed fringes to extract for each channel several time-specific 
values for amplitude, bit-stream alignment (BSA) delay, fringe phase and phase- 
delay rate. In the process, fringe phase is corrected by subtracting the phase 
obtained from calibration tones in order to remove unwanted instrumental effects. 

The phase and BSA delay extracted from all channels can be used to obtain a 
final measured delay by means of bandwidth synthesis. This process first com- 
putes ambiguous delays by means of pairwise combinations of fringe phase from the 
different frequency channels. Ambiguities are then removed by a multistep proc- 
ess that begins with the unambiguous BSA delay and iteratively removes ambiguities 
in the BWS delays for ever-wlder channel-pair separations until the widest separa- 
16n (and most precise delay) is reached. 
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Fig. 1. Simplified illustration of antenna instrumentation 









For each observation of a given source (1-10 minute duration), the output 
observables of the BLKO systen are usually one composite value per channel for 
each of the following: fringe amplitude, fringe phase, phase-delay rate and BSA 

delay. Further, one value of BWS delay is obtained for ea'-h channel pair. From 
these observables, one can extract a best value of delay and of delay rate for 
subsequent processing. Currently, the "best" delay is usually the BWS delay 
from the most widely spaced channel pair. When two RF bands (typically 3.6 cm 
and 13 cm) are recorded, the final value for delay (delay rate) is obtained by 
combining the best values for delay (delay rate) from the two bands in a way 
that removes the effect of charged particles. The resulting delay should be 
equal to the sum of the geometric delay and the clock offset, while the delay 
rate is the time derivative of that delay. 

A block diagram of the data reduction steps found in the current version of 
tbe BLKO system is presented in Appendix G. This report omits later data reduc- 
tion steps that simultaneously fit the delay and delay-rate observables from many 
observations to obtain estimates for geophysical, astrometric and clock 
parameters. 
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SECTION III 
CLOCK MODELING 


Before proceeding to cross-correlation analysis, it will help to define a 
few time-keeping terms. First, define 

t = true (universal) time 

= true oscillator frequency (relative to true time) 

w = nominal oscillator frequency (e.g., 5 MHz) 
c 

Two "clocks" will be defined to assist in the analysis. The time of the 
first clock, the "recorder clock", is registered by the bits on the tape. This 
bit time is given in terms of true time by 


S = ^ ^ "b 


where t, is the error in the recorder clock. 

D 


( 1 ) 


The second clock, the "(>-cal clock", will be established, at least concep- 
tually, near the point where the 5-MHz signal from the station enters the phase 
calibrator Instrumentation, as illustrated in Figure 1. The phase at this refer- 
ence point is given in terms of true frequency by 


= 0) (t - t ) (2) 

c t o 

where t is the (unknown) true time at which 4>(. = 0. As implied by Eq. (2), we 
will assume that the oscillator is perfectly stable so that the only Imperfections 
in the station clock are offsets in epoch and rate. If phase (t»j, is measured at 
the reference point, (Ji-cal time can be defined by 


t 

c 




(3) 


where we assume t^ is zero when (|>(. is zero. This is possible since a particular 
zero crossing can be defined to be = 0. If the nominal frequency is related to 
the true frequency by 


0) = (I) + Au) 

c t c 


then 0-cal time is related to true time by 


(A) 


t 

c 


t + T 

C 


(5) 
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where the clock error is given by 


T M 

C 


A(i) 



t ) - t 
o o 


( 6 ) 


Usually the difference between and is 


small, less than a part in 10^. 


It is conceptually and instrumen tally convenient to define a clock in terms 
of a 1-pps signal that drives a counter to register the number of received pulses. 
To establish the l-pps signal, a particular set of position-going zero crossings 
(one every 5 x 10^ crossings for 5 MHz) must be selected. (Equivalently in 
Eq. (2), a definition of (|ic - 0 must be made.) This arbitrary selection for the 
4>-cal clock leads to an ambiguity which must be resolved in BLKO clock synch 
measurements, as discussed in Appendix G. Clock synchronization with the BLKO 
system will be measured relative to the ^-cal clocks. 
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SECTION IV 


THE RECORDED NATURAL SOURCE SIGNAL 


This section derives an expression for the analog signal by tracking the 
Incoming signal from Its free-space form through the Instrumentation to Its final 
form at baseband. In this and subsequent sections, only a nonrelatlvlstlc theory 
will be considered. 

' " 

Id^VLBI measurements, two antennas simultaneously receive the signal pro- 
duced by a very distant, very compact natural source of random radio emissions. 

In this section, we will assume for simplicity that the natural source Is a point 
source and leave It to Appendix A to generali'ze the analysis to an extended 
source. (That generalization greatly complicates the analysis but changes the 
final cross-correlation function only by Inserting an additional factor In the 
amplitude and an additional term In phase.) Due to the ^reat distance of the 
natural source, the Incoming signal can be modeled by plane waves. Typical 
natural radio sources used In JPL VLBI work emit noise with a wide smooth fre- 
quency distribution. We can model one member of the random ensemble of noise 
waves as a superposition of plane waves In the form 


E(x,t) 



A((a>) exp[l((A)t - k • x)] do) 


( 7 ) 


where E Is the electric field at true time t and point x, k « u>k/c Is the apparent 
wave vector of the plane wave and gives the apparent direction of propagation, 
and A(ai) Is the (random) Fourier amplitude at frequency u). All quantities are 
measured with respect to a geocentric frame with axes defined by true equatorial 
coordinates of date. The adjective "apparent" refers to the aberration effect 
that will be observed In a geocentric frame. For simplicity, the wave is 
assumed to be linearly polarized. We will also assume at this point that the 
wave propagates in a vacuum and will not include propagation media effects. 
(Troposphere effects can be included as a simple additive group delay and would 
add little to the analysis. Space-charge effects are treated in Section XIV.) 

Relative to the geocentric frame, the electric field at a given station 
becomes 


Ej(t) = E(xj(t),tj 





( 8 ) 


where xj(t) is the location of station j as a function of true time. Station 
position Is defined to be the Intersection of axes of the antenna. (We will 
Include below an Instrumental delay Tq to account for the fact that the natural 
source signal Is Injected at the horn rather than at the intersection of axes.) 
Phase Is expressed In terms of true time since that time will be common to both 
antennas In a nonrelatlvlstlc analysis. 
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After Injection, the signal passes through several filters and Is 
heterodyned to baseband under the combined effect of various mixing signals^ 

The phase effects of all Intervening components can be grouped Into three cate- 
gories. First, the overall effect of the mixing signals can be described by one 
total mixing signal, o>hjl^j ♦hj* where Is the total mixing frequency and 

represents both a constant phase offset and variations In mixing phase that 
are nonlinear In time. Second, all group delays after the Injection point. 
Including effective group delays through filters, can be represented by one total 
delay tjj. Third, all Instrumental phase shifts (except and those of the 
form wTjj) can be represented by one total shift We also Include a term 

for delays before the tone Injection point. This delay Is defined to be the 
difference In path length for the radio wave In Its actual propagation to the 
tone Injection point and Its theoretical propagation to the Intersection of 
axes. This delay Is necessitated by the definition of station location and 
Is theoretically Incorporated by delaying the Incoming plane wave. The formula- 
tion with Tg assumes that dispersive phase shifts above the Injection point are 
negligible. If they are not, then additional calibrations will be required to 
remove them. 

Figure 1 schematically Illustrates the modifications to signal phase as the 
signal progresses through the Instrumentation from injection to recorder. The 
block diagram has been simplified to show only one element of each type with 
token phase shifts and group delays. Nevertheless, the following analysis is quite 
general if the variables ace properly defined to represent lumped effects of all 
the elements in an actual system. 

When the lumped instrumental effects are included, the doppler-shif ted analog 
signal at baseband for a given BWS channel at station j at bit time tbj can be 
written in the form 

^ 00 

Vj(t^j) = / A(w) Gj(yj) e^'^j do) + c.c. + (9) 

where 

y^ = o)(l - k • Xj/c) ^10) 

♦j - - "aj> - ^ ' “h) 'j ■ *hj - 1i ’ij ■ 

In this expression, Gj is the amplitude of the system bandpass, nj is additive 
noise, ya Is the doppler-shlf ted (station) frequency corresponding to u> and 
c.c. denotes complex conjugate. The time argument tj is the true time correspond- 
ing to bit time tjjj (see Eq. 1). 

In the BLKO system, the signal is sampled at a 4-Mbit/s rate and recorded in 
a bilevel mode. That is, one bit is recorded for each sample point and that bit 
Is 1 if the signal is positive and 0 if the signal is negative. In subsequent 
analysis, we will ignore this bit code and act as though +1 or -1 were directly 
recorded on the tape. 
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Before proceeding to cross-correlation analysis » we will reformulate the 
expression for the analog signal ao all components are normalized In the mea- 
sured units of noise temperature. For station J, let the normalized analog sig- 
nal corresponding to V be given by 







( 12 ) 


where z, q, n are the total signal, natural source signal (first two terms In 
Equation 9) and additive noise (Iasi term in Eq. 9), respectively. We will 
assume these quantities are all nomullzed so that 


^z^^ ■ ^q^> “ 1 (13) 

The quantities T^, T., Tg are the noise temperatures for the total signal, the 
natural source and additive noise respectively. It is assumed that the spectral 
power for each noise source is fairly constant across the 2-MHz passband so that 
one average noise temperature can describe the whole passband. (The noise power 
per passband is actually kTgW, but the factor kW drops out in the ratios.) As 
we shall see, the normalized form in Eq. (12) is more convenient for cross- 
correlation analysis. 
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SECTION V 

THE CROSS-CORRELATION FUNCTION 


The BLKO correlator (the CIT/JPL Mark II correlator) aligns the bits on the 
two tapes, multiplies them bitwise and fringe-stops (counter-rotates) the result- 
ing fringes. This section develops a detailed model for the cross-correlation 
function, the expectation value of the product* of the aligned signals. We will 
explicitly treat only one frequency channel, with the understanding that the cor- 
relator also processes the other channels in their time sequence. 

The expectation value of the bitstream product, the cross-correlation function 
(fast fringes), is defined by 

r(t„.o ^ - tb) (14) 

where Tj^^, are model delays used to offset the two bitstreams (see Appen- 
dix B). We will use a tilde to denote a digital or quantized version of a quan- 
tity as for the two-level value (*1) for Vj and for the lag (tqj^), which is quan- 
tized to the nearest bit. The term £b explicitly represents arbitrary lag off- 
sets (£ - -7 to +8) about the central lag of the correlator, by offsetting bit- 
stream i in units of the bit interval b. For the BLKO system, b is equal to 
(4 Mbit/s)“l or 250 nsec. The argument tQ is the common Greenwich time defined 
by the bits on the two tapes. That is, the correlator implicitly enforces the 
assumption that 


‘■bi 


bj 


= t. 


(15) 


and uses t^ as the time tag. 

The van Vleck relation (Ref. 2) gives the cross-correlation function in 
Eq. (14) in terms of the analog signals: 




<Vi(t 


bi - ’.1 - » ■ >> 




- '.i>> 


<V^> 


(16) 


*In the BLKO correlator, one bit stream is first multiplied by a stopping 
sinusoid and the resulting prc.''uct is multiplied by the other bitstream, as 
explained in Appendix B. Since multiplication is associative, we will reverse 
the actual order for clarity. 
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where small correlated amplitude (r << 1} has been assumed. (In typical VLSI 
measurements, one finds r < 0.2.) With the normalised representation of the 
preceding section, the expectation of the normalized voltage product in Eq. (16) 
becomes 




(17) 


:v‘Xv‘> 




(18) 


since the additive noise nj is uncorrelated between stations. The arguments 
t),j are shorthand for the sums of terms that are the arguments in Eq. (16). Thus 
the problem is reduced to finding the correlation between the normalized natural 
source components (qj^) for the two stations, with the overall amplitude scaled 
by a multiplicative system temperature factor. Note that this factor times 
2/it gives the amplitude of the cross-correlation function for perfectly corre- 
lated natural source components (i.e., <qiqj> ■ 1.0). 

We will now present on the basis of a plausibility discussion the final form 
for the cross-correlation function. A complete derivation is given in Appendix A, 
including the effects of source extent and the random nature of the source signal. 


The natural source signal q is given by the first two terms in Eq. (9), 
normalized to have unit power. Thus the product q^ qj is given by a double 
Integral over frequency where the integrand Includes the sum of two products of 
the form A(u))A*((n' ) and A((ii)A(u'). Under the expectation-value operation, only 
th product <A(u))A*(ui' )> has a non-zero value and this occurs as a delta function 
aL 0 > 0 )’. (For a stationary random signal such as that from a natural source, 
the amplitude at one frequency is uncorrelated with the amplitude at another 
frequency.) Thus, the double integral over frequency collapses to a single 
integral over fr'i^quency with one term in the Integrand. If we let Sp((i>) denote 
the source power spectrum arising from <A(u))A*(ai)> . the cross correlation func- 
tion becomes 



(19) 
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where 


r = V 


Sp(u)) G^(yj)du)^ Sp(u) Gj(yj)du. 


( 20 ) 


( 21 ) 


= ♦h “ “hj "mj ■" “hi 


+ k 


[*j(tj - 


W “ ■ "«1 ■ 


+ w(T. +T.+T -T.+T,“ £b) +♦-+♦_ 

b I a mi mj IB 


( 22 ) 


in which 


* ♦«> • <"hi'i * ♦hi* 


(23) 


and Ttj, T 3 , and are differences between the corresponding station quanti* 
ties in the order j minus i. 


The terms in the cross-correlation function in Eq. (19) can be explained as 
follows. Normalization of the signal in terms of system temperature generates the 
factor rg,, as suggested by the derivation of Eq. (18). Included in rg, is the 
fringe visibility Vf, which is equal to the correlated flux divided by the total 
flux. This factor and a brightness transform phase shift three 

modifications to the final cross-correlation function when the analysis is gener- 
alized to Include extended sources (see Appendix A). The third modification is to 
change the point-source wave direction £ to ^ 3 , a wave direction defined relative 
to a somewhat arbitrary reference point within the extended source. The factor 
of 2/i> is the amplitude loss due to two-level sampling (see Eq. 16). The factor 
of 1/2 is a consequence of complex representation of the fringes and gives the 
proper definition for the power spectrum Sp (i.e., by cancelling the factor of 2 
gained when the complex conjugate in Eq. (l9) is added to obtain a real number 
for the cross-correlation function). 

J 

In the denominator, the quantity Du is the product of the total source power 
from the two stations and arises from the normalization of the natural source 
components at each station (that is, from the requirement <q^^> ■ ^qj^^ " 1 ). 

Note that, as one would expect, the total power at each station is equal to the 
Integrated power across the passband where the power at each frequency is given 
by the spectral power of the natural source multiplied by the square of the band- 
pass amplitude. 

Within the Integral in Eq. (19), the various terms can be explained as 
follows. Since a component of the natural source signal at one frequency is 
uncorrelated with a component at any other frequency, the correlated signal can 
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be obtained by multiplying each frequency component at one station by the 
correaponding component at the other station and summing over all frequencies. 

By corresponding frequency components, we mean the components with the saaw 
frequency in geocentric coordinates. Due to the earth's rotation, a given geo- 
centric frequency will be doppler shifted when measure by station instrumenta- 
tion. Those shifts are accounted for by the doppler-shif ted arguments yj^ in the 
bandpass functions G|(. 

Since interferometers are differencing devices, the phase in Eqs. (19) 
and (22) l.s the difference in the phase for the two stations, each with the 
proper oi'^'s-jt in time. All of the phase terms in Eq. (22) can be derived in this 
manner except for three which Involve special considerations. First, the bright- 
ness transform phase is a consequence of source structure as shown in Appen- 
dix A and is included without proof. Second, in the instrumental delay terms 
yj Tij, we have neglected small doppler effects of the order 

(k • x./c) T-. ^ lO”^ * lO”* sec • 10 sec (24) 

a J IJ 


which can shift the phase by 0.01 cycle at X-band. In practice these snull 
doppler terms can be removed by applying a correction based on instrumental group 
delay, if necessary. Thirl, the quantity Tg arises from the difference tj^ - tj 
which, through the use of Eq. (1) plus (15), can be shown to equal 


- t, 


‘bj 


- T 


bi 


(25) 


where is the total error in bit time at station k. 

Introduction of t 1, an ^.nterim delay approximately e^ual to the usual 
geometric delay, simplifies fringe phase We define Tg by 


« 


UT 


g 





i 


'.1 - “>>] 


(26) 


This delay is discussed in Section XII. Note that tI depends weakly on i, so 
that X-band phase changes by less than 0.005 cycle when one lag is traversed. 

In the following analysis, we neglect this dependence and evaluate Xg at the 
lag with peak fringe amplitude. Since this lag dependence can be accurately 
modeled, it can be easily accounted for during processing, if necessary. When 
Eq. (26) is substituted into Eq. (22), fringe phase becomes 



(27) 


The phase rate if can reach about 20 kHz at X-band. Therefore, before the data 
can be compressed by summing over bits, the fringe phase aust be counter-rotated 
to a very. low rate (<0.1 Hz). The next section discusses that procedure. 
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SECTION VI 


FRINGE COUNTER-ROTATION 


In nost observations of natural sources* the systen noise is much greater 
than the correlated amplitude so that the SNR for a single bit product (Eq. 19) 
is very small (0.0002 - 0.2). Thus many bits of data ('^'10^ - 10^) imist be col- 
lected to raise the SNR to a usable level. However, due to the high phase rate 
that the fast fringes typically possess (as high as 20 Ulz at X-band), the data 
volume cannot be compressed by directly summing the fringes over time (bits). 
This section explains fringe stopping, a process that counter-rotates the 
fringes (i.e., digitally heterodynes the fringes to a low frequency) and thereby 
allows sums over an interval of time (correlation interval) that compress an 
observation containing many bits ('^10° to 10^) into a relatively small number 
(100 - 1000) of points (correlation coefficients). In addition to permitting a 
massive data compression, counter-rotation can, with modification, partially 
remove the effects of model-delay quantization, as we shall see. 

For fringe stopping to be successful, the model phase and model delay used 
by the correlator mist accurately reproduce the time dependence of the actual 
fringes. Specifically, to avoid significant losses in amplitude in the output 
fringes, the geometric model must be accurate to a part in 10^ or better. The 
model phase can then be constructed with sufficient accuracy to insure that only 
a small fraction of a cycle of the stopped fringes will be traversed in a corre- 
lation interval (e.g., 20 kHz x 10~^ x 1 sec ■ 0.02 cycle). Thus, amplitude 
loss due to incoherent addition can be avoided. We will assume the correlator 
model delay is sufficiently accurate. 

The model phase* used by the correlator for counterrotation is given by 



(28) 


wh^ere is the best estimate of and u>* is the best estimate of the centroid 
(ui) of the bandpass product. The time t^ is a reference time adopted ar the cor- 
relator and is usually the start of the observation. The last three terms involv- 
ing quantized delays are designed to r^iove the average quantization effect from 
the overall stopped phase (at effective frequency 


To counter-rotate, the voltage product discussed in Section V is aailtiplied 
at each bit by the appropriate model phasor and sumsed over N^ bits within a 
correlation interval, producing one complex value for the fringes for that 
interval : 


“«“o> ^ T 52 - *•> V‘c - 


(29) 


^For simplicity, we will neglect a constant- frequency term that can be incluefed in 
model i^se. 
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where E denotes a sum over the bits In the Interval and t^ Is now the common bit 
time at the middle of that interval. (As explained above through Eq. IS, tQ is the 
Greenwich time defined by both bitstreams.) The sum Interval (correlation inter- 
val), which ranges from 0.2 seconds to several seconds and contains at least 
8 X 10^ bits, is equal to the total dwell time for the given BWS channel within 
a multiplexing cycle. By combining Eq. (14) with the expectation value of 
Eq. (29), one finds that the expectation value for the complex stopped fringes 
is given by 


t I 


When Eq. (19) is substituted for r(t,)l), two terms result: a sum-frequency term 

and a difference-frequency term. For sufficiently high fringe rates 
(^f ^ 100 Hz), the sum-frequency term essentially averages to zero in the sum 
over time leaving only the difference term: 


<u,(t,)> 


-1-ijV f‘s e 0 


d(i) 


where the stopped phase is given theoretically by 
A0f = Of - V 

= A(|i. + u)At + - (w - (D. .) ib + (o) - u)' ) (t - t ) + (fr 

h I nl m m o 

for which we have made the definitions 


At = t' + T. + + T - T 

g b I a m 


(31) 


(32) 

(33) 


(34) 


‘♦h = ♦h - <“'hj 


- (t 


bl “ 's* 


(35) 
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mi 


- T 


mj 


(36) 


Two terms of the form (‘*^hj“*‘'hj order of 10”^ cycle or less 

have been omitted. We would emphasize that Is the model delay used by the 
correlator at the central lag (1 = 0). 


The fringe-stopping analysis in this section Ignores the fact that the BLKO 
correlator actually uses a trllevel quantized model for the stopping sinusoids, 
as explained in Appendix B. We will assume that the use of such a model changes 
the amplitude of the stopped fringes by a multiplicative factor but otherwise 
will not significantly change Eq. (31). That multiplicative factor is the 
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amplitude of the first harmonic of the quantized stopping sinusoid and is equal 
to 1.176 for the quantization scheme used in the BLKO correlator (see Appendix B). 
One can see how this amplitude factor enters by expanding the quantized sinusoid 
in terms of its harmonics and showing that only the first harmonic beats the 
cross-correlation signal to a low frequency. That is, the sum over correlation 
interval filters out all terms except this low frequency term. Thus, except for 
the coefficient of the first harmonic, the resulting fringes are assumed to be 
the same form as those generated by the pure stopping function of this section. 
Exceptions can occur in special cases due to interactions between various fre- 
quency components in the signal and stopping function, but those cases will not 
be treated In this report. 

The theoretical expression for the stopped fringes can be rewritten by 
interchanging the sum and IntegriJ in Eq. (31) to obtain 



(37) 


where 


= A(|), + u)At + (j) + (J> 

I h 1 o 


(38) 


F (w - Ci)' ) = 

q 



0)') (x 

Dl 


T ) 
m 


(39) 


The factor a^ has been inserted to account for quantization of the stopping 
sinusoid, as discussed above. The sum over time has been applied only to the 
"quant ized-delay phasor" (from fifth term in Eq. 33) since the time variation 
of that factor will dominate when an accurate model delay is used. That is, 
the rate of the residual phase (At|/f) can be made very small (<0.01 Hz) so that 
the phase Atlif will change by only a small amount (<0.01 cycle) over a sum 
interval (''-1.0 second). On the other hand, the phase (w-io') (x^-Xm) changes 
by 0.25 cycle over a delay quantization update interval, which can be as 
small as 0.1 second. As explained in Appendix D, the quantization filter Fq 
can easily be calculated on the basis of delay quantization history. Exceptions 
to the above assumption of small residual rate can occur in "initialization" runs 
involving a poorly known baseline or source position. However, this will occur 
only for the first pass through the correlator since the results of that pass 
can be used to improve the a priori for all subsequent passes involving that 
baseline or source. 


In the BLKO correlator, correlation sums are made over time intervals that 
are often large compared to the delay quantization update interval (bit-jump 
interval) . An alternate correlator design is to form intermediate sums for very 
short intervals over which the delay quantization error changes onj.y by a small 
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fraction of a bit. For each of these shorter intervals « the correlator irould per- 
form a Fourier transform from the lag domain to the frequency domain. After the 
transformation, the quantization error appears purely as a phase effect that can 
then be removed by appropriately counterrotating fringe phase. This approach 
produces a slightly better SNR (by 0 to 3.5%) than the BLKO system because the 
correlated signal is added coherently for each frequency across the 2-HHz 
passband. However, it also increases the complexity and cost of the correlator. 

By accepting a small loss in SNR, the BLKO system can onploy a simpler corre- 
lator and let postcorrelation software handle the quantization problem with 
the quantization filter Fq. Even though the quantization problem Itself does 
not necessarily require it, the postcorrelation software (l.e., the phase 
tracking program) does eventually transform to the frequency domain for other 
reasons, as discussed in the next section. 
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SECTION VII 


TRANSFORMING TO THE FREQUENCY DOMAIN 


After the stopped fringes leave the correlator, they are passed to the 
phase-tracking program (PHASOR) . In the present version of PHASOR, the stopped 
fringes are Fourier transformed from the lag domain to the frequency domain. 

The reasons for this transformation are as follows. In the lag domain, the 
phase is implicitly referenced to the true centroid frequency, which depends 
on bandpass shape. Since this reference frequency enters in the calculation 
of BWS delay, uncertainties in bandpass shape can introduce errors in BUS 
delay. In the frequency domain, an exact reference frequency can be forced 
on the data so that output phase will not depend on bandpass shape. Thus, 
frequency-domain analysis eliminates a possible source of error and avoids 
the operational burden of precisely determining bandpass shape, at least for 
phase measurements. However, for the most precise measurements of correlated 
amplitude , it may be necessary to determine bandpass shape accurately. 

One disadvantage of the frequency domain approach is that, in practical 
applications, a Fourier transform can introduce phase errors. Those errors arise 
because in practice only a small number of lags can be summed in the transform, 
leading to poor resolution (» bandwidth/4) in the frequency domain. However, 
accuracy studies show that the errors introduced are not as serious as one might 
guess. When complex fringes based on typical bandpass shapes were simulated (see 
Appendix E) , it was found that a nine-lag transform reproduced the actual overall 
input phase to better than O'. 1 millicycle (with the peak amplitude placed at the 
central lag to within 0.5 bit). For a nearly rectangular bandpass, the "odd” 
frequencies (n/2N, n odd) more accurately reproduced the input phase than did 
the "even" frequencies. The transformed shape deviated from the input shape 
by as much as 9% for the 9-lag transform. These amplitude deviations should be 
no problem since the transformed shape can be modeled if the actual bandpass 
shape is known. An alternate approach is to determine directly the transformed 
shape by measuring for a strong source of known strength the amplitude of the 
cross-correlated signal in each frequency bln of the frequency-domain fringes. 

Two other considerations arise in the choice between domains — phase callbra 
tion and fractional-bit-shift correction. Both of these corrections can be 
carried out in either domain, although they are more easily applied in the fre- 
quency domain. Thus, primarily because of the decoupling of phase from amplitude 
frequency domain analysis appears to be the preferred approach in phase- 
tracking and that approach has been followed. 

The remainder of this section transforms the fringes from the lag domain 
(Eq. 39) to the frequency domain. For reference. Appendix C reformulates the lag 
domain fringes to cast them in a form used in the past. 

Before transforming to the frequency domain, one must detect the fringes 
and determine which lag offset produces the largest fringe amplitude. This is 
accomplished by performing an FFT over time on the fringes for each lag offset 
within a lag range believed to contain fringe power. The output of these FFT's 
can be viewed as a function of lag and frequency (where the frequency variable 
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is fringe frequency associated with time dependence rather than observing 
frequency). The amplitudes of the FFT points are searched over lag and fringe 
frequency to determine which lag and which frequency bin has the largest 
amplitude. The maximum amplitude is checked to verify that it exceeds the SNR 
threshold. Only fringes for lags near this central lag are used in subsequent 
processing. 

The stopped fringes for each correlation Interval (Eq. 37) are transformed 
to the frequency domain by a discrete Fourier transform over a subset of the 
16 lags passed by the correlator: 





ioi/fb 

k 


e 


(40) 


where b is the Interval, Iq is the central lag with the largest amplitude 
and «.* is summed over lags centered about Hq (Nil = N2 - Ni + 1). In the 
present BLKO system, a 9-lag transform is used since the phase errors associated 
with a nine-lag transform are apparently the best for values of Nji less than 12 
(see Appendix E). For Nj^ = 9, we would have = -4 and N£ = 4 . The frequen- 
cies are baseband frequencies spaced at equal Intervals across the system 
passband. To obtain an adequate sampling rate in the frequency domain while 
maintaining statistical Independence of the noise on the transformed points 
h(d)j^), the frequencies are given the spacing 


k+1 
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As observed above, the "odd" frequencies give the best phase accuracy so we will 
use 
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(42) 
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(43) 


where W is the nominal channel bandwidth of 2 MHz (half the sample rate). For 
N|^ ■ 9 only the fringes associated with the first four values of k (i.e., k ■ 1 to 
4) have significant amplitude. The frequency values for these points are W/9, 
3W/9, 5W/9, 7W/9. The point at W (k ■ 5) has very small amplitude for the BLKO 
passband. 


Some of the comments found earlier in this section can be clarified by calcu- 
lating the expectation value for the transformed fringes: 
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where where the plus (minus) sign Is used for USB (LSB). The sign 

selection depe^s on the sign of u-uhi* which is positive for upper sideband (USB) 
and negative for lower sideband (LSB). By substituting the expression for <U|i> 
found in Eq. (37) and interchanging the sum and integral, one obtains 


<h(w|^)> 
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where 
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One can show 
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(46) 

(47) 

(48) 


where 6(x) is the delta or Impulse function. Thus, for an infinite number of 
lags we would obtain precisely the cross-correlation spectrum: 




There are practical limits to the number of lags that can be transformed. 

For example, as more lags are used, the power is spread over more points in fre- 
quency, making data reduction more costly. Further, the cost of the correlator 
increases as the number of lags Increases. On the other hand, transformation 
errors decrease as the number of lags increases. A reasonable compromise for 
the BLKO system appears to be 9 lags. For this value, the power is spread over a 
manageable number of points (4) in the frequency domain while transformation 
errors are fairly small, as discussed in Appendix E. 

In the following analysis, we neglect the small transformation errors and 
use Eq. (51) to describe the phase in the frequency domain. If necessary, the 
transformation errors can be corrected to a large extent through modeling. 

It is informative to note that, for a flat power spectrum Sp and for 
Identical bandpass shapes at the two stations, the bandpass ratio in Eq. (49) 
can be converted through the use of Eqs. (21) and (50) to the form 


S 

_E 




G^.k> 


(53) 


where the bar denotes an average across the passband. We have assumed that the 
passband amplitude is nearly zero at W (2 MHz for the BLKO system) and that 
doppler shifting is negligible (Af^j/W % 0.01 for the BLKO system). For a 
rectangular bandpass this ratio is 1.0. 

With the fringes in the frequency domain, the uncertain bandpass amplitude 
appears as a multiplicative factor. The fact that this form (Eq. 49) decouples 
phase from amplitude can be demonstrated heurlstlcally by computing phase at 
each frequency by means of the inverse tangent of the ratio of real and imaginary 
parts. (This method is not actually used for phase computation, as seen in 
Section X.) In such a computation, the amplitude would cancel in the ratio. 

As mentioned above, another advantage of the frequency domain is that phase cor- 
rections can be directly and simply applied to the fringes by means of a phase 
counter-rotation to each frequency component h(u)^). That technique is used in 
the phase calibration procedure discussed in the next section. 


SECTION VIll 


STOPPED TONE PHASE 


This section digresses to treat the subject of phase calibration by means 
of the tone generator and derives an expression for stopped tone phase. Appen- 
dix F derives expressions for tone amplitude and tone SNR. 

In VLB! measurements, the observed phase Is corrupted by unknown and 
unstable shifts due to Instrumentation. Such phase effects can degrade the 
accuracy of geophysical measurements and complicate measurements of clock 
synchronization. However, most of these Instrumental effects can be removed 
through the use of a phase calibrator. The commonly used approach, pioneered 
by A. E. E. Rogers (Ref. 3), Is to Inject at a point near the front of the 
Instrumentation a calibration signal consisting of a set of tones, equally 
spaced In frequency and derived from the station frequency standard. Properly 
used, this technique can both correct for system phase Instabilities and allow 
absolute calibration of Interferometer phase so that clock synchronization Is 
possible. 

The calibrator signal Is Injected at a low power level and Is Imbedded In 
the radio source data. The phases of the Injected tones for the current tone 
generator are expected to be stable at the lO-psec level and can be calibrated 
absolutely to the nanosecond level. Each channel bandpass will contain more 
than one tone, but the exact number In a given application will depend on system 
characteristics and specific accuracy goals. At the correlator, each tone Is 
counter-rotated to nearly zero frequency by digitally mixing each single-station 
bitstream with a frequency that closely approximates the baseband frequency 
for that tone. After the resulting tone phase has been extracted and differenced 
between stations. It Is subtracted from Interferometer (fringe) phase. Since the 
tones are subjected to the same Instrumental effects as the natural source signal 
Instrumental effects after the Injection point cancel In the final difference. 
(This report will not consider difficulties encountered In calibration of Instru- 
mental effects before the Injection point.) After instrumental calibration, the 
measured delay Is reduced to the sum of a geometric delay and a clock synchroniza 
tlon offset, as we shall see In the following sections. (For simplicity, propa- 
gation media terms are Ignored until Section XIII.) 

In this section, tone phase Is separated Into the various effects that enter 
the calibrator signal In transit from the clock through the Instrumentation and 
then through data reduction. Since the analysis will be based on an "ideal" 
model for the system, extra effort might be required to account for deviations 
from ideal behavior. It Is beyond the scope of this report to assess all compli- 
cations. 

The model for stopped tone phase is developed In the following sequence: 

(a) from the ^-cal clock to the Injection point, (b) from the injection point to 
the recorder and (c) through the data reduction procedure. The simplified block 
diagram of the Instrumentation shown In Fig. 1 will be useful in tracing the sig- 
nal through the Instrumentation. 
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First, for step (a), the phase of the signal at the clock reference point 
((^-cal clock) discussed in Section III and illustrated in Figure 1 is given in 
terms of true time and nominal frequency by 


as Indicated by F,qj. (3) and (5). We have chosen to represent calibration phase 
in this way since (a) the clock error Tj. is explicitly shown, (b) true time t is 
common to both stations, and (c) tone frequencies are calculated at the correla- 
tor on the basis of the nominal frequency u>c> The signal will experience a 
cable delay so that its phase entering the tone generator will be given by 

4>c • (i)^ <t + - T^) (55) 

The "ideal" tone generator detects he zero crossings of the 5-MHz input signal 
and converts the positive-going zero crossing into rectangular pulses as indi- 
cated schematically in Fig. 2. The nominal width (fp) of the output pulses for 
the current tone generator is about 20 psec. Repetition rates slower than 5 MHz 
are obtained by blanking pulses. For example, by passing every 10th pulse, the 
pulse repetition rate would be 500 kHz. The passed repetition rate will be 
denoted by oip. As Indicated below, the tone frequencies are the harmonics of 
the fundamental repetition rate cop. 

In practice, there will be deviations from the ideal waveforms shown in 
Fig. 2, deviations such as delays, phase shifts and an amplitude distortion 
arising in both the tone generator and connecting elements. For example, a zero 
crossing is not precisely at the center of a rectangular pulse and the pulse 
shapes are not perfect rectangles. This report will not consider such deviations 
h”.t most of these effects can be formally Included by making a complex Fourier 
expansion rather than the cosine expansion presented below. (It should be 
pointed out that a constant deviation from rectangularity will cause clock syn- 
chronization measurements to be biased by a constant error that Is less than the 
pulse width but will caur;e no error In geophysical/astrometric measurements.) 

The output of the tone generator can be decomposed into its harmonics, which 
will have frequencies given by niUp (see Appendix F). The phase of the n*^^ calib- 
ration harmonic (tone) at the injection point will be given by 


= U) (t + T 

n n c 




(56) 


The small group delays from the input of the tone generator to the injection 
point have not been explicitly represented since they (or any other group-delay 
between (^-cal clock and Injection point) can be included in t^. 

We are now prepared to model effects on tone phase Introduced in transit 
from the injection point to che recorder and then through the tone-stopping 
process. After injection, the tones, along with the natural source signal, 
pass through various filters and mixers until thev reach baseband within a 
2-MHz bandwidth. As in the case of the natural source signal (Eq. 11), the 
effects of these instrumental stages can be represented by lumped variables so 
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thet the phase of the tone recorded at bit time t|,j for station J 
becomes 




U) (t. + T . 
n J cj 
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( 57 ) 


where the instrumental phase effects are given by 


.t, + . + <** ) 
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Note that the delay does not appear since the phase-calibration signal is 
not subjected to antenna delays. Again, tj is the true time that corresponds 
to bit tiiw tjjj at station j. 

The recorded tones are processed at the same time the tapes are cross- 
correlated. The correlator operator will use a best estimate for the baseband 
frequency of each calibrator tone to separately "fringe-stop” each tone in a bit- 
stream (see Appendix F). That is, in effect, the correlator will subtract from 
the phase of the n^^' tone the phase 


♦nj ■ <“n - “hj> <'bj - 
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(59) 

(60) 


where is the best (nominal) estimate for the heterodyne frequency, tuj is 
bit time'^at station J, and tg is the reference time adopted by the correlator. 

For phase calibration to work properly, it is necessary to use the same values 
for (i)f,j, and (jj^j, and tg in both the tone model phase (Eq. 59) and the cross- 
correlation model phase (Eq. 28). Based on Eqs. (57), (58) and (60), the theo- 
retical expression for the resulting stopped-tone phase at station j becomes 

A(|i (uj ) - (Ji - (61) 


which is equal to 
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where the "stopped" heretodyne phase is given by 
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(63) 


This puts the stopped-tone phase in the form we vant, expressed in terms of con- 
tributing effects due to instrumentation and data reduction. 
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This analysis has neglected some subtle but important problems that can 
seriously degrade the stopped tone phase. These problems arise from unwanted 
Intermodulation of frequency components in the calibrated signal with frequency 
components in the stopping '.unction, where the resultant, possibly aliased 
frequency of the Intermodulation product happens to be close to the frequency of 
the stopped tone. With careful design, these problems can be substantially 
avoided. 

Several methods exist for extracting tone phase from the stopped tones. The 
method used in the present version of the BLi^ system is discussed in Appendix G. 
In the following section, we will assume tone phase has been successfully 
extracted. 

Another fact ignored to this point is that measured tone phase, like 
interferometer phase, has integer-cycle ambiguities. Ambiguity resolution for 
tones falls into two categories — relative and absolute. By relative ambiguity 
r^'solution, we mean the removal of ambiguities in a relative sense between tones 
within a channel passband. If such ambiguities are present, the variation of 
phase with frequency within a passband cannot be determined and correct inter- 
polations in frequency cannot be made. For example, if one wishes to obtain an 
absolute estimate of bitstream-alignment (BSA) delay after phase calibration, 
one must either explicitly or implicitly determine the correct effective phase- 
frequency slope for the tones in the passband. To remove relative ambiguities, 
one must have an a priori estimate of the slope of phase vs. frequency, which 
means one must have estimates for terms multiplying in Eq. (62). Those terms 
are the correlator reference time ts, the calibrator cable delay* (Ty), the in- 
strumental delay from the tone injection point to the recorder (t^), and the 
epoch offset (it, - t ) between the recorder clock and the ^-cal clock. The max- 
imum allowed a priori error on these terms depends on tone spacing. For example, 
if the closest tones are separated by 500 kHz, then the collective a priori error 
(3 o) on these terms must be less than one microsecond. Ambiguities on the other, 
more widely spaced tones can be removed by an Iterative procedure akin to the 
bandwidth synthesis method. An example of relative ambiguity removal is given in 
Appendix G. 

After relative ambiguities have been removed, there will remain in tone 
phase an overall (absolute) ambiguity which is the same for all frequency points 
in the passband. Thus, the overall calibrated interferometer phase obtained from 
each passband will not be corrected for absolute integer cycles due to instrumen- 
tation. As we shall see, it is not necessary to consider absolute instrumental 
ambiguities. Accurate measurement and application of only the fractional part of 
overall Cone phase is sufficient for the ultimate removal of absolute ambiguities 
in corrects J BWS delays and corrected interferometer phase, provided the BWS 
channels are properly spaced. 


More exactly stated, T^ Includes all delays from the clock reference point to 
the tone injection point. 
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SECTION IX 
PHASE CALIUUITION 


The next step in data reduction is to correct the fringe phase from each 
BWS channel with the phase derived from the tones injected by the phase calibra- 
tor. We will asstnne that there are several tones spaced across each passband 
and that integer-cycle anbiguities in tone phase can be rmsoved in a relative 
sense between tones (see Appendix G) so that calibration phase at any frequency 
in a passband (i.e., at the frequencies in Section VII) can be obtained through 
interpolation. We will neglect errors in this process such as interpolation errors 
that might occur as a result of phase ripple across a passband. 

The calibration phase at frequency (D|( within a passband is given by the dif- 
ference in tone phase obtained for the two stations: 

Based on Eq. (62), this phase is theoretically equal to I 

♦t<v • ^i> ^i ^♦h 

where i|,» tui tj, ^^h differences between stations in the order 

J minus i. For simplicity we have assumed that all dispersive phase shifts in 
In Eq. (il) occur after the tone injection point so that they appear iden- 
tically in tone and fringe phase. If significant dispersive phase shifts occur 
before the injecticn point, separate calibrations will have to be applied. 

For each frequency component h(iD^) in the transformed fringes (Eq. 49), we 
will assume tone phase is subtracted from fringe phase by sieans of a counter- 
rotation (multiplying by e“^^t) of the fringes at each fringe point (correlation 
interval). Other methods can be used to apply phase calibration at various 
stages in processing, but all of those approaches will not be discussed here. 

Fringe-point counter-rotation is chosen for the analysis in the text because it 
purifies the phase at an early point, making the fringes easy to interpret and 
use at an early stage. However, phase calibration in the BUCO system is currently 
applied to the final overall value obtained for each observable from an observa- 
tion. That appr.iach changes the order of the operations described here and leads i 

to essentially the same final value for all observables. See Appendix G for a j 

description of that procedure. 

When the fringes in Eq. (49) have been counter-rotated with the tone phase 
in Eq. (63), one obtains 


( 66 ) 
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where the corrected phase is given by 


‘♦c ■ “k <’i * ’c * ■ 'u ■ ♦«‘“k> 


(47) 


We have dropped the lag-shift term in Eq. (49) since it is known and can 

also be counterrotated out of fringe phase along with tone phase. Removal of 
this term references the model delay to the central lag >*f the processor, thereby 
insuring that the residual phase from all BWS channels is based on the same model 
delay. As discussed in Section X, the current version of PHASOR does not correct 
for the lag-shift rerm until phase tracking of each BWS channel is completed. 

We have removed it here to simplify the analysis. 

At this point, it is instructive to make a simple interpretation of 
interferometer phase given by Eq. (67). If separate measurements were made to 
determine the cable delay and antenna delay t^, these quantities could be 
removed from phase. The model delay Xg| is known exactly and could also be 
removed. (In practice, these quantities are actually removed at a later stage 
in processing, as discussed in Section XIII.) The interferometer phase would 
then be left in a clean and useful form: w|^(xg ). This is the phase that 

would result from measurements made at RF by ideal instrumentation (l.e., no 
cable delays, phase shifts, instabilities, etc.). Such ideal instrumentation 
can he visualized as a "point RF recorder" placed at the intersection of axes 
with no instrumental components to complicate the phase that is received and 
directly recorded at RF at each station. Under such ideal assumptions, one 
can easily derive the interferometer phase (ig -f Tc) by means of a simple 
subtraction of the RF phases from the two stations after offsets to account for 
geometric delay and clock offset. 

In the calibration process described above, a subtle point has been 
neglected. Only the fractional part (in cycles) of both interferometer phase 
and tone phase can be determined when phase is extracted after counterrotation. 
Thus, when tone phase is subtracted from interferometer phase, the subtraction 
actually denotes removal of the fractional part of instrumental phase. We can 
ignore integer-cycle errors at this point since we realize that there are 
Integer— cycle ambiguities due to delay modeling that must be resolved later in 
the BWS process. The BWS process requires as input precise val'tk x for only the 
fractional part of interferometer phase, free of unwanted instrumental phase 
(fractional part). The integer-cycle part (due to and x^.) is then determined 
in a recursive process described in Appendix H. The output is a delay measure- 
ment free of postinjection Instrumental effects. It is a peculiarity of this 
technique that, even though BWS delay and interferometer phase can be ultimately 
calibrated absolutely with regard to instrumental effects, instrumental phase, 
including integer cycles, is never absol 'ely measured. 
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SECTION X 


PHASE TRACKING 


At this stage, we are prepared to extract the observables of amplitude, 
phase, phase-delay rate, and BSA delay from the stopped fringes. More than one 
method Is available for extracting observables from Che fringes. The two most 
familiar approaches are phase-tracking and coherent Integration across channels. 
The coherent Integration approach (Ref. 4) Is based on relatively simple, fast 
algorithms and has the Important advantage of better SNR for signal detection. 
However, for this approach to work, the fringes must be made coherent across all 
summed channels by means of phase calibration. Since the phase-tracking approach 
does not require such 'Coherence, we chose that approach in the early development 
of the BLKO system wh> a phas“ calibration was not available. Further advantages 
of the phase-track app^ ach Include (a) versatility In modeling and weighting 
fringes, (b) easy access to intermediate observables for system trouble- 
shooting, (c) computation of RMS fringe residuals, and (d) simple adaptability 
to single-channel spacecraft applications with poorly known nonlinear phase. 

In the phase-tracking technique, the fringes are analyzed (by a program 
called PHASOR) to extract all of their Information content, including phase as 
a func‘-ion of time during an observation. To obtain the time dependence, each 
observation Is divided Into equal sections, whose length falls In the range of 
a few seconds to a few minutes. As discussed below, the transformed fringes are 
fit within each section with a parametrized function to determine four quanti- 
ties: amplitude, bit-stream-alignment (BSA) delay, phase-delay rate, and phase. 

The connection of phase between sections Is discussed below. Phase tracking 
results in a further compression of the data in that each observation will be 
characterized by four parameters for each of 1-30 sections Instead of 200-1000 
complex fringe values. 


As suggested by Eq. (67) and demonstrated below, a suitable mathematical 
model for the fringes within a section is given by 


T(o.j^,tg) = A^ g(u,^)Fq e^^ 


( 68 ) 


where the model phase is given by 


= ( 0 )^ - 
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(69) 


and 



+ 0 ), 


hi 


(70) 


The parameters Ag, Xg, x^jig and (|>g are the four solve-for parameters mentioned 
above. The complex filter F^ accounts for the delay quantization and is derived 
in Appendix D. The reference frequency Wg and reference time tg will be placed 
at the center of the bandpass and center of the fitting Interval, respectively, 
in order to minimize correlations and errors for phase, delay and delay rate. 
With this selection, the solve-for phase 4ig will be the phase at frequency 
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and time As implied by the fitting equation, it is assumed that the fringe 

phase can be adequately modeled as a linear function of time and as a linear 
function of frequency over the fitting Interval. For some applications, the 
amplitude function g(ui^), which becomes a real four-point array called BAMPL in 
PHASOR, can be calculated with adequate accuracy from the equation 


g(»;> 


2 

— a, = — 

^ 1 r2 


( 71 ) 


where G is the measured bandpass amplitude (assumed to be the same for both sta- 
tions as in Eq. 53). For a rectangular bandpass, we obtain g(ii>^) ° 2ai/v. For 
more accurate amplitude measurements, a more complicated relation than Eq. (71) 
can be used as discussed in Appendix E. 

Since the model is not linear in all parameters, a priori estimates of the 
four solve-for parameters must be obtained before a least-squares fit can be 
performed. Three of the parameters (Ag, ig, x^g) are obtained from the FFT's 
over time (one for each lag) that were used to detect the central lag (see Sec- 
tion VII). FFT Integration times can range, in theory, from a few seconds to an 
observation length, depending on the application. The output of the FFT's over 
time can be viewed as a function of lag and fringe frequency (ideally with a 
sin x/x dependence in both directions) that peaks for unique values of those 
quantities. The determination of the location of this peak in lag-frequency space 
gives a direct estimate of the BSA delay (xg) and residual fringe frequency (and 
thus delay rate, T^g). The peak amplitude gives an estimate of fringe amplitude 

^s- 


Once Ag, Xg and x^g have been estimated, one can make an initial estimate 
of (|)g by means of a least-squares search (with the frequency-domain fringes) that 
minimizes the differences in the actual fringes and the model fringes as a func- 
tion of 4ig. In this process, the model fringes are computed from Eq. (68) using 
as fixed values the estimates of Ag, Xg and x^g explained above. In the search, 

<|)g is incremented in small steps (''^0.1 cycle) from 0 to 1 cycle. (This procedure 
will reveal only the fractional part of phase. Procedures for dealing with the 
integer part will be discussed below.) For each phase Increment, the square of 
the difference between the actual and model fringes is summed over time points in 
the given section. The resulting values for the (''^10) sums are then searched for 
the minimum that reveals which phase value was the best estimate. (In practice, 
an Interpolation between steps is carried out to obtain a more accurate estimate.) 

After all four a priori estimates have been determined, the fringes are 
subjected to a standard least-squares fit that simultaneously estimates all 
four parameters for the section. Theoretical expressions for the resulting 
parameter estimates can be formally demonstrated as follows. Comparison of 
Eq. (66) and Eq. (68) indicates that the solve-for amplitude will be given by 
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Further, the solution will require the phases to be equal so that we will have 
from Eqs. (67) and (69) 
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where 2nTr represents the Integer-cycle ambiguity mentioned above. Expanding the 
right-hand side of this equation in terms of frequency and time about and tg 
and equating like terms on each side, one finds that the solve-for parameters 
will be given by 


(74) 


(75) 


(76) 


where 



At, =t’+t +t -t -t 
t g c a u m 


(77) 


All of these terms are evaluated at time and frequency uz- The delay Tq is 
the model delay at the central lag of the correlator, as indicated in Section VI. 

Thus, all three of these solve-for parameters depend on Atf. Note that 
(ignoring the 2nTr ambiguity and structure effects in ^B) = 4>s/uz = Axf, while 

delay rate is approximately equal to the time derivative of that quantity. This 
relationship between ts and (fis is a fundamental requirement for single-observation 
bandwidth synthesis (BWS), a process that is used to account for effects of 
integer-cycle ambiguities in The removal of ambiguities between frequency 

channels by single-observation BUS is discussed in Appendix H. For each BUS 
channel, ambiguities are removed between sections within an observation by means 
of another procedure to be discussed below. 

Examples of fringes in the frequency domain are given for an intercontinental 
baseline in Fig. 3 and for a short baseline in Fig. 4. The fringes in Fig. 3 are 
the frequency-domain counterparts of the lag-domain fringes presented in Appen- 
dix C. The solve-for values of amplitude, frequency, phase and 6SA delay are 
given at the top of the figure, where the errors are the covariance errors based 
on RMS fringe residuals computed for the fit interval. (Note that amplitude has 
been multiplied by 10^ for convenience.) Except for amplitude, the solve-for 
values in Fig. 3 are in agreement with the values for the same parameters esti- 
mated in the lag domain in Appendix C. Amplitude is different since different 
normalizations are applied in the two domains. That is, the solve-for amplitude 
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Long-baseline example of stopped fringes for the BLKO system 


6 








REAL COMPONENT 


i 






In the lag domain is equal to the amplitude of the fringes li: Eq. (37) while. 

In the frequency domain. It Is given to first approximation by Eq. (72). The 
current ratio for these normalizations Is about 2.0. For these particular fits, 
the amplitude normalization matrix, BAMPL, has been assigned the values (1.02, 
1.14, 1.14, 1.14) for the four frequency bins. These values appear to match the 
data fairly well and are not very different from the theoretical BAMPL values 
calculated In Appendix E on the basis of nominal bandpass shape. However, the 
differences between the two BAMPL arrays Illustrate the level of difficulty In 
obtaining one bandpass shape that Is suitable for all data. (For very accurate 
amplitude measurements. It may be necessary to use a strong source of known 
strength to generate calibration fringes from which both the relative magnitudes 
of the BAMPL components and the overall BAMPL normalization can be determined 
for each channel.) 

By Inspecting the plot, one can readily see that the solve-for values are 
in accord with values easily estimated by Inspection. The frequency can be 
estimated from the period of the fringes. The phase Is referenced In time to 
the center of the interval (arrow) and in frequency to a nominal centroid 
frequency (0.9 MHz). A value for phase can be estimated by Interpolating to the 
reference frequency at the reference time. BSA delay can be roughly estimated 
by extracting the shift in phase (see dotted arrows In Fig. 3) between the two 
outer frequency bins In the plot and dividing this shift by the frequency dif- 
ference (2W/3 = 4/3 MHz) between the bins. 

The short baseline fringes in Fig. 4 are provided to illustrate the dis- 
continuity in the fringes introduced by a particular 1-bit jump in delay, as 
modeled through Fq (see Appendix D) in Eq. (68). As mentioned earlier, such 
discontinuities are also present for intercontinental baselines but occur 
so frequently that a number of jumps are averaged over in a single fringe point. 

Up to this point, we have considered parameter estimation for only one 
section In an observation. In most cases, an observation is divided into several 
sections of equal length (5 - 120 seconds), each of which Is separately fit for 
As* Tg, Tjtg, and 4>g. An observation can be divided into sections so that there 
is less chance for the time dependence of the solve-for parameters (particularly 
nonlinear phase rate) to violate the assumptions underlying the fringe model and 
so that the Mme dependence can be observed. However, the section must be long 
enough to provide enough "SNR” for a successful fit. 

Since fringe fits yield only the fractional part of fringe phase, addi- 
tional effort is required to properly connect phase between sections. The phase 
for the first section Is forced to fall In the range 0-1 cycle. (It is this 
arbitrary assignment that gives rise to the absolute ambiguity in overall phase 
from the observation.) A preliminary estimate of the phase of the second section 
is made by adding to the least-squares phase of the first section the phase change 
predicted by multiplying the estimated fringe frequency (from the FFT or from 
fringe fitting) by the time Interval between section centers. This preliminary 
estimate is then compared with the more precise but ambiguous phase produced by 
the fit to the second section. The two are forced Into agreement by Increment- 
ing or decrementing the fit phase by Integer cycles. Clearly, the phase 
extrapolation to the second section must be in error by no more than ±0.5 cycle 
for this procedure to work. Once the ambiguity in the second section has been 
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removed, its phase and the fringe frequency are used in the same fashion to 
resolve the relative ambiguity in the third section, and so on, until the last 
section is completed. 

Figure 5 displays the phase values (lower curve) associated with the fringes 
in Fig. 3. However, to obtain those phase values, the fringes in Fig. 3 were 
rerun through PHASOR with the section length reduced from 60 seconds to 
12 seconds. This shorter section length was chosen to better reveal the time 
behavior of phase. As explained in the figure, a phase correction equal to 
(<*>k ~ <^hi) ^o^ has been applied in order to shift the measured phase, (^s((*>z), from 
the peak lag (Hq) to the reference lag of the correlator (1*0). (See Eqs. 49 
and 66 and the comments after Eq. 67.) Note that the phase rates of the two 
plotted channels are about the same, as one would expect, and that the nonlinear 
phase excursions of the two channels are highly correlated, as one would expect 
for many error sources (e.g., ionospheric fluctuations and station oscillator 
instability). 

In summary, the output from this stage in processing consists of values 
for amplitude, BSA delay, phase-delay rate and phase, computed for each frequency 
channel as a function of time (section) across an observation. A formal error 
for each parameter value is computed by means of a standard covariance analysis 
based on the RMS residuals found in the fringe fits. 
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SECTION XI 

DELAY MODEL IMPROVEMENT 


Although it Is adequate for cross-correlation, the delay model used by the 
correlator Is fairly crude and Is Inadequate for subsequent processing. After 
phase tracking, PHASOR replaces the Inaccurate correlator delay model with a 
more accurate model by appropriately adding back to the three solve-for param- 
eters given theoretically by Eqs. (74-76) the correlator model (which is 
numerically passed) and then subtracting a new more accurate model. As indi- 
cated by Eq. (26), the geometric component of the new model should accurately 
reproduce the "geometric" delay that has been set up In fringe phase by time- 
offsetting in the correlator. The new more accurate model delay is given by 
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where we have used the central lag (1 q) for the lag offset I under the assumption 
that such an approximation is sufficiently accurate or that lag-dependent cor- 
rections have been made (see Section V). The station positions ^jjjj(t) are highly 
accurate models for the motions of the stations, while k„, is an accurate model 
for the source direction. In the most demanding applications, one must make a 
distinction between the bit time (tbj) and true time (tj) when calculating the 
above delay, since the geometric delay will be parametrized in terms of true time. 
This can be accomplished by inputting a clock error (xbj) for each station and 
using it to compute true time from bit time (see Eq. 1). The resulting values of 
true time are then used in model calculations as suggested in Eq. (78). The con- 
fusing presence of a different true time for each station is a consequence of 
correlator design, as discussed in Section V. 

Model replacement is carried out for each of the three aforementioned solve- 
for parameters from each section in an observation. For each section, the time 
arguments in the model delay are determined from the section reference time (i.e.. 
tbj = tbj^ = tz = section reference time so that t^ ■ tg - and tj ® t^ - 
Tbj ) • 


In addition to the geometric delay, one can also replace at this step the 
model for troposphere, the clock or other delay components if desired. 

In subsequent sections, we will still refer to Eqs. (74-76) as the theoreti- 
cal forms for the observables, but with the understanding that the model delay 
has been replaced with an improved value. The resulting change in the observa- 
bles can be accounted for by changing Eq. (77) to the form 
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wliere Tq is the new model delay. BSA delay, BWS delay, and phase delay (except 
for cycle ambiguities and structure effects) are now all theoretically equal to 
this expression for delay but are known at different levels of precision. 

Even though model improvement is always applied in actual data reduction, 
we will suppress it in our example observation so that data flow can be continu* 
ously traced from figure to figure. 
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Fig. 5. Example of stopped fringe phase 
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SECTION XII 


OBSERVABLE VORHATION 


After the four solve-for quantities (unplltude, B8A delay, phase-delay 
rate, phase) have been extracted from the fringes for each section across an 
observation, they must be analyzed to obtain the final observables. As men- 
tioned above, the behavior of the parameters as a function of section (tlin) 

Is often informative when there are equlfment problems. Ionosphere signatures, 
and the like. However, to reduce final data volume. It Is usually deslr^le to 
pass on to subsequent software one composite value per observation for each 
parameter. Composite values for wnplltude and BSA delay are obtained in PHASOR 
by averaging the values for all sections in an observation. On the other hand, 
all of the delay rate values from Individual sections are discarded since a 
better estimate of delay rate can be derived from the section phases as follows. 

As illustrated in Fig. 5, the phase values at the section centers of a BUS 
channel are fit by least-squares with a linear function of time to generate a 
composite phase-delay rate (slope) and a composite phase (constant). The compo- 
site phase from this fit becomes the final estimate for the RF phase for that 
channel. In order to minimize correlations, the reference time for the coaiposlte 
phase Is defined In the fit to be the center of the observation. After the, fit, 
phase residuals are computed to reveal nonlinear trends. In addition, a covariance 
analysis based on the formal errors (obtained from fringe fitting) In section 
phase Is performed to obtain the formal errors in the composite phase and phase- 
delay rate. 

As mentioned above, the usefulness of the phase observable is greatly 
reduced by an absolute ambiguity that currently remains unresolved in the BLKO 
system. On the other hand, the composite phase-delay rate from the fit is 
absolutely known (except for occasional and resolvable aliasing) and becomes the 
composite observable for the given channel. For the final observables of phase- 
delay rate and BSA delay, one has the option of either using the composite value 
from one channel or averaging over all channels. 

The calculation of the composite BWS delay observable is smoewhat aiore 
complicated. First, BWS delays are computed for each section of a channel pair 
by combining the section phases according to the formula 
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where the phases (^a.<^b) frequencies (u>a,(i>b) come from channels (a, b). From 
Eq. 76, one finds the theoretical value for BUS delay is given by 
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where the partial of structure phase approximates the associated finite difference 
(assuming that Aw/w <<< 1) and where the last term represents the integer-cycle 
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aabigulty. To roveal tlM depoodoaee acroos an obaervation, a aaparato MS dalay 
vnlua la conpucad for aach aaction in an obaarvatlont an lllustmcod in Fig. 6 for 
tha phaaa values in Fig. S. Int^tr-cycla wbiguicias are then reaoved in a 
relative sense by comparing delays from adjacent sections. (Residual delay 
typically changes very slowly across an observation.) The section delays from 
each channel pair are then fit by least^squares with a linear function of time 
(aoction) to determine a composite BUS delay rate (slope) and a composite BW8 
delay (constant) for that pair. In the fit, the reference time is defined to be 
the middle of the observation. The BWS delay rates for the sepsrate channel pairs 
should be equal, within errors, to one another and (ignoring the effects of space 
charge) to the phase-delay rates obtained above for the separate channels. Since 
the BWS delay rates are less precise, they serve mainly as a consistency check. 

(In cases with very poor a priori delay m^els, BUS delay rates can be used to 
reiBove aliasing in phase-delay rates.) It should be pointed out that the rates 
for the particular case in Figs. 5 and 6 should differ by about IS sHs due to an 
artificial frequency-independent offset inserted in model phase during fringe 
stopping. If that offset, which was included for fringe plotting purposes, is 
removed from the rate in Fig. S, the rates from the ttn> figures are in fair 
agreosent. 

Along the final composite values, we have an ambiguous BUS delay from each 
channel pair and an unambiguous BSA delay. It is clear from Bqs. (74) and (81b) 
that all of these delays are given theoretically by Eq. (77) plus the source 
structure partial, except for integer-cycle ambiguities In the BUS delays. That 
makes it possible to eliminate the BWS ambiguities by single-observation band- 
width synthesis, a bootstrap technique that begins with the unasd>iguous (but 
least precise) BSA delay and successively removes ambiguities in the (increasingly 
precise) BWS delays for ever-wider chaimel separations. The details of this 
procedure are given in Appendix H and are illustrated by cos^leting the reduction 
of the long-baseline observation presented in preceding sections. When phase 
calibration is available, a slight improvement in precision can be obtained by 
using an alternate method to obtain BWS delay. In that method, the channel phase 
values, corrected for relative ambiguities, are fit with a linear function of ob- 
serving frequency to obtain the BWS delay (slope). After extraction by either meth- 
od, the BWS delay is given theoretically by Eq. (81b) without the ambiguity term. 

In suamary, the processing of a single observation produces one composite 
value for each of four observables: amplitude, phase-delay rate, BWS delay and 

SF phase. Ambiguities are absolutely removed in the BWS delay but are not cur- 
rently removed in RF phase. Therefore, RF phase is not used in subsequent process- 
ing for geophysical inforsmtion although it probably will be in the future. 
Presently, the composite BWS delay from the most widely spaced channel pair is 
generally taken as the final delay observable. 
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SECTION Xlll 

CALIBRATION OP PROPAGATION HEDU DELAYS 


For simplicity, we have neglected in the preceding analyals the delays 
due to propagation media. In the case of the troposphere, the propagation delay 
enters simply as a group delay (i.e., adds another term to Mf in Bq. 80). Con- 
sequently, once the troposphere has been calibrated, its effect can be removed 
from delay or delay rate by a simple subtraction. On the other hand, the treat- 
ment of delays due to charged particles is a little more complicated. 

The effect of charged particles on fringe phase (in radians) can be repre- 
sented by 


K N 

♦e<“> • ( 82 ) 

where u is the observing frequency (rad/sec), iL is a constant, and N^ is the 
difference in the integrated electron content along the two raypaths of the 
interferometer. The group delay corresponding to this phase shift is given by 
the usual expression: 


Tg(w) 



K N 


c e 


2 

li) 


(83) 


It can readily be shown that, neglecting source structure effects, charge- 
particle effects will show up in the interferometer observables (Eqs. 74-76) in 
the form 


At, + T (u) ) 
f e z 


(84) 

Af j i^(u.^)/u.^ 


(85) 

ATf ^ ♦g(o)^)/u)^ 

(neglecting 2nir) 

(86) 


where is the reference (centroid) frequency for the channel under considera- 
tion. In some ex,>erlments, the final BSA delay and delay rate are obtained by 
respectively averaging Tg and f^g over all channels in the band. In these cases, 
the frequency Ug in Eqs. 84 and 85 should be changed to the effective frequency 
given by 


-2 


ui 


eff 



(87) 


where N is the number of channels, 
c 
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For the two-channel BWS delay (Eq. 78), the charged-particle delays enters 
In the form 


'‘"bus = ‘'t 


*^(u) ) - 
e a e b 

CU - 0), 

a b 


( 88 ) 


a At. - 


K N 
c e 


f 0) 0). 

a b 


(89) 


Thus the effective frequency for the charge-particle delay in two-channel At 
is given by 


BWS 


= ^ 


BWS 


a b 


(90) 


We are now in a position to discuss S/X calibration of charged-particle 
effects. Suppose that we have measured the interferometer observables in 
Eqs. 84, 85, 86 and 89 at both S-band and X-band. Note that ts, 

^BWS equal to a group delay term plus a charged-particle term proportional 

to (d~ 2. One can easily show that the o)“2 term for any one of these observable 
types can be removed by means of a linear combination of the S-band and X-band 
observables for that observable-type: 


0 = c 0 + c 0 

SX S S XX 


where 0 and 0 are the observables at S- and X-band and 
s X 


c 

s 


2 // 2 
= (i) / ( 0) 

s s 


2 


01 


X 


) 


(91) 


(92) 


c = - 01 ^) (93) 

X X X S 

The effective frequencies uig and o>x are computed for each observable type as dis- 
cussed above. With this choice of the constants, the calibrated observable 0 ^^ 
will be independent of charged-particle effects (i.e., only the first term in 
Eqs. 84, 85, 86 or 89 will remain). 


The observables at S- and X-band can also be used to extract the electron 
content (or its rate). The constants used in the linear combination for this 
purpose are 



2 2 ,. 2 

01 01 / (ui 

S X X 



( 94 ) 
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This second linear combination will give 


O' - -N, 
sx 3 

from 3SA delay t 

8 

(95) 

- N 

e 

• 

from delay rate t. 

98 

(96) 

“«e 

from phase delay ^ /u (neglecting 2 nir) 

S 2 

(97) 

» -N 

from BWS delay 

(98) 


We should note that, if residual observables are used in the S/X combination, 
as in the BLRO system, the model delay implicitly applied in At£ must be the same 
for the S-band and X-band observables. In presenting the expression for the 
coefficients in Eqs. 92, 93, and 94, we have assumed combinations of one obser- 
vable type: BWS delay with BWS delay, phase delay with phase delay or delay rate 

with delay rate. If BWS delay from one band is combined with phase delay from 
the other, then the coefficients will change somewhat. Those expressions can be 
easily derived and will not be presented here. 

After S/X calibration, the residual observables are theoretically free of 
charged particle effects except for dispersive errors in measured delay. One 
"dispersive" effect in dual-band calibration, the source structure terms 
neglected in Eqs. 84-86, is discussed in Reference [5]. 
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SECTION XIV 
MOIffiL RESTORATION 


For computational and interpretive convenience, residual delays (i.e., 

Eqs. 74, 75, 76 and 81b) rather than total delays have been utilised in the 
reduction up to this point. However, since the desired final observable is the 
total delay with the model delay restored, one must conq>ute the model delay that 
must be added to the residual BHS delay. The geometric conq>onent of this model 
delay is chosen so that the geometric component of the resulting total delay is 
equal to the conventional geometric delay. The conventional geometric delay is 
defined as the difference in time of arrival of a wavefront at the intersection 
of axes of the two antennas, where one station is taken as reference. As 
derived in Appendix I, the model geometric delay resulting from this definition 
is given by 


r (t ) 
mg 1 






c(l + k 


m 




(99) 


where station i has been chosen as the reference station and where the 
denominator accounts for the motion of station J during wave transit (transit 
time correction). The source position ^ and station-position ^i(t) in this 
expression must be exactly equal to those used in the model removed (Eq. 78). 
The time argument ti, which becomes the time tag for the final observables, is 
set equal to true time at the middle of the observation. (Since residual delay 
and delay rate essentially remain the same over time intervals that typically 
separates true time from bit time, one has the freedom to make this choice for 
time tag). 


It might seem strange that the restored model (Eq. 99) is not equal to the 
model removed (Eq. 78). The models differ in the offsets in the time arguments 
of the station positions (compare Eq. 78 with Eq. 12 in Appendix I). In the 
model removed, both stations must be given argument offsets in order to match the 
phase produced by the correlator when both bitstreams are offset by model delays. 
The modeled argument offsets for the two stations must be precisely equal to 
those correlator model delays? In the model restored, the conventional defini- 
tion of geometric delay requires that only station j be offset in time where 
chat offset is equal to the geometric delay (see Appendix I). The reason one has 
some freedom in selecting the restored mojlel delay is that the delay effects 
due to such argument offsets are small (igij Tg/c < 2 x 10~^ * 20 msec > 40 nsec). 
This means that they can be accurately modeled and, therefore, accurately 
removed or inserted. He can then choose to restore the model that makes the 
final total delay conform to the conventional definition of geometric delay. 


The total model delay to 
Tq except the geometric model 


be restored is equal to Tgg plus all other terms in 
Tmg (see Section XI) and is most easily expressed as 


mA 


T + T* - T* 

mg m mg 


( 100 ) 


After this model is added back and after measured values for the "cable delay" 
Tu and the "antenna delay" Ta are removed, the resulting delay observables from 
Equations 76 and 81b (with Atf given by Equation 80) become 
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. . 2nu 

• T + — 4- 

♦ “z “* 


( 101 ) 


for phase delay and 


^Bws * ^ air 


( 102 ) 


for BWS delay (with ambiguity removed) where 


T = T + T 

g c 


(103) 


The total geometric delay, which entered the preceding equations in the form 


T = t' - t' + T 

g g mg mg 


(104) 


now has the conventiowl definition given in Appendix 1 and is defined relative to 
a reference position k^. Thus, except for structure effects and an integer-cycle 
ambiguity, the BUS delay and phase delay are equal to the same delay: a geometric 

delay plus a clock synchronization term. When these observables are obtained 
through S/X calibration, the structure terms are evaluated at X-band. Strictly 
speaking, the structure effect is a geometric effect and can be absorbed in the 
geometric delay by placing the reference position at the effective position [5]. 
However, a single reference position cannot be used to make the structure effect 
disappear for both and Tgyg since the effective positions for the two observa- 
bles usually are not the same. 

For computational convenience, phase-delay rate, which as a measured residual 
has the form given in Equation 75 (but is obtained from residual phase slope) , is 
passed to subsequent software not as an observed rate but as two delays whose 
difference will represent the delay rate. The measured residual phase-delay rate 
is converted to two total delay values by the relations 

r,„(t + 6t/2) = T .(t + 6t/2) + At. 6t/2 (105) 

9 K mA 9 

T,„(t - 6t/2) = T .(t - 6t/2) - A^ . 6t/2 (106) 

tnA 9 


where 6t ~ 1 to 4 seconds. In subsequent processing, these delays are converted 
to a delay-rate-like observable by the relation 


6t “ «t 

This reconstructured rate observable is essentially equal to 

. • ^B 

T . = T + — 

9 U)^ 


(107) 


( 108 ) 


which is the time derivative of the phase-delay observable. 
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With regard to structure effects, it can be shown (5), that, for each obser- 
vation, one effective position can be found that, when used as the reference 
positiM ka* structure terms in both BUS delay and phase-delay rate 

equal to zero. Thus, those observables can be represented by Tg -f Tc and fg ■<- 
wtere Tg is the geometric delay relative to the effective position for those two 
observables. 
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SECTION XV 
SUMH^Y 


For each of several baselines (^10) » the BLKO syston can cross-correlate 
and compress approximately a gigabit of recorded data for a natural source to 
four final observables. After phase-calibration and dual-band charged-particle 
calibration, the four final observables are fringe amplitude (Equation 72), 
phase-delay (Equation 101), phase-delay rate (Equation 107), and BWS delay 
(Equation 102). Multiple baseline observations of fringe amplitude, perhaps 
coupled with triangle-closure values for the other observables, can be used in 
subsequent analysis to study source structure. If the delay and delay-rate 
observables are measured for many sources, they can be subjected to a multi- 
parameter fit to obtain astrometric and geophysical parameters. 

The Block II system now under development will be superior to the BLKO system 
in many ways. The Improvements include (a) a lower fringe-detection threshold, 
as well as simpler and faster operation, through coherent integration across fre- 
quency channels (a BLK II option that requires phase-calibration) , (b) more fre- 
quency channels (28), (c) a higher total bit rate through continuous recording 
of each frequency channel, (d) better phase calibration by using four tones per 
channel and by frequent ''real-time'' counter-rotation of fringe phase (at every 
correlation Interval), (e) reduced intermodulation and aliasing in stopped tone 
phase by quantizing the tone-stopping sinusoids to 128 amplitude levels rather 
than 3, (f) improved capability for spectral-line work (224 lag integrations), 

(g) longer integration times for a given data-storage-array size by means of 
data compression through digital filtering, and (h) an improved fractional-bit- 
shift correction through frequent ('v> every 10 msec) counter-rotation of fringe 
phase. The Block II system is scheduled for completion in 1982. 
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APPENDIX A 

A CCMPLETE DERIVATI(»I OF THE CROSS-C(HtRELATI(H( FUNCTION 


This appendix presents a complete derivation of the cross-correlation 
function, taking Into account the statistics and the spatial distribution of 
the missions of the natural source. 

The radio noise generated by a very distant extended natural source can be 
expressed as a superposition of plane waves in the form 


E(x,t) 



exp[loj(t - k-x/c)] doi dfi + c.c. 


(Al) 


where E(x,t) is the electric field at position x and true time t. For stoplicity, 
the electric field is assumed to be linearly polarized. The function A(k,u) is 
the Fourier amplitude at frequency u for the wave received from direction Ic, while 
c.c. denotes complex conjugate. As in Section IV, all q\}antities are measured 
with respect to a geocentric frame. The wave direction k can be expressed as a 
function of two parameters and, in the case of right ascension and declination, 
becomes 


k = -(cos6 cosa, cos6 sina, sin6) (A2) 

where a, 6 are the apparent right ascension and declination relative to true 
equatorial coordinates of date. The quantity dD represents a differential solid 
angle such as cos6dad6 in the case of right ascension and declination. In the 
following steps, we will use the wave vector; k = wk/c. 

The electric field detected at antenna j is given by 


E^(t) = E(Xj(t),t) 



gi[(i)t-^*Xj (t) ] 


db) dD + c.c. 


(A3) 


where *j(t) is the position of antenna j as a function of true time. (See Sect- 
tion IV for a definition of antenna location.) This signal will be received by 
the antenna and passed through various filters and mixers. The modeling of the 
composite effect of instrumentation is discussed in Section IV. Following that 
derivation, the signal component due to the natural source can be written in the 
form 




J j A(k,b.) Gj(y^) dD + c.c. (A4) 
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where the doppler-shifted passband argument Is given by 


ITj ■ 0)(1 - k-Xj/c) 


(A5) 


and where Instrumental phase terms are given by 

= “hj 'j * * J'j ’ij * 

As in the text, tj is the true time that corresponds to bit time tj^j. 


(A6) 


We will first compute the expectation of the product Of two analog signals 
as described by Eq.'(A4) and then compute the normalization. The expectation 
value of the voltage product is given by 

I 


^V*^(t ) v"(t 



^A(k,o))A (k' ,o)')^Gj^(y^)Gj(yj)e^'^ do) dfi dw' dfi* 


Iflf 


+ I I i I <A(k,(o) A(k' ,u)' )>G^(y^)Gj (ype^'^^c dw dfi du' dJJ' 


+ c.c. 


(A7) 


where 


= w (1 - k*Xj/c) 
yj = w' (1 - k'*Xj/c) 


\l/ = wtj^ - k*x^(t^) - oj'tj + k'*Xj(tj) - a^(u),t^) + Oj(u)',tj) 


(A8) 

(A9) 

(AlO) 


and \l>f. Is a similar expression that will not be needed. The quantities t)^ are 
the offset values for true time corresponding to tjji^ as in Eqs. (16) and (17) 
in the text (i.e., tj,j * + A implies tj “ tj + A). 


We will, assume that the natural source is completely incoherent, which 


means 



(All) 
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where Sp(kyu>) Is the power spectrum for direction k and 6(e) is a Dirac delta 
function.^ That is, anissions from different areas of the source are uncorre- 
lated. Furthermore, the noise emitted by a given area is stationary and there- 
fore possesses uncorrelated frequency components. Since the signal is real, we 
know 

A(k,u)) - A*(k,-w) (A12) 

This relation and Eq. (All) give 

<A(k,u>)A(k',o)')> « i Sp(k,o)) 6(k-k’) 6(u+u)’) (A13) 

* 0 for u) and u* > 0 (A14) 


Under these assumptions, Eq. (A7) becomes 

[r Sp(k,u)) Gj(yj) exp(li|»j)da) dD + c.c. 


where 


" -V 

u) (1 - k*x./c) 


0) (1 - k'Xj/c) 


<l'i = w + k*B^ - a^(u),t^) + Oj(o),tj) 

The A(k,u)) A(k',u)) term has dropped out since u> - -u>' is not covered in the 
region of the integration. We have also defined a "retarded baseline:" 


(A15) 

(A16) 

(A17) 

(A18) 


” *i^*^i^ (A19) 

Let two particular parameters (S,y) define the direction vector. Suppose 
the brightness distribution is very narrow about some reference direction k^ 


k 

a 




(A20) 


For two particular direction parameters (S, y), the delta function denotes 
6(6 - B*) 6(y - y'). Furthermore, we will require B and y to satisfy the 
relation dBdy ■ dW. 
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If the brightness distribution is sufficiently compact, we can approximate Che 
wave direction by 


k ■ k. + 96 


<« - «.> * I7 


(V - y.) 


(A21) 


a 


where the partials* are evaluated at the reference point $ , y • Ne then obtain 

d A 

<vj(tbi) V^(tjjj)> “ J I R(u,v,«i)) G^(yj) Gj(yj) e*p(i<l»j) + c.c. 

•'a 


where 


y^ ■ 10 (1 - k^‘Xj/c) 


yj - (i) (1 - k^*Xj/c) 


<l»f *= u) • a^(o),t^) + aj(u),tj) 


(A22) 

(A23) 

(A24) 

(A25) 


and where 


R(u,v,a)) 



Sp(e.Y.w) 


^-2..1(n(e-e^)+y(y-y^))jjj, 


for which 


u = - 


3k 

3B 


. 3k 
‘ 3 y 

X = 2tic/ui 



(A26) 


(A27) 

(A28) 

(A29) 


*Structure coordinates are usually defined as 6 ■ (a > 0 |.) cos dc and y " ^ 
where (oc, 6^) is the origin of structure coordinates relative to r.a. and 
declination. With this definition, the 6 partial and y partial are in the 
direction of increasing right ascension and declination, respectively. 
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with this definition* u and v are the componente of the baeellne vector 
projected on the plane of the sky. Two approxlnatlons have been made In 
Eq. (A22). First, the weak ic dependence In the bandpass functions has been 
neglected so that the yi*yj values have been evaluated at Second, the limits 
of the (6 ,y) Integration have been extended to Infinity under the assumption 
that the brightness distribution is very compact and terminates the Integration. 
The function R(u,v,w) Is referred to as the brightness transform. A more 
detailed treatment of the brightness transform and of structure effects Is given 
in Ref. 5. 

The fringe visibility Is the normalized modulus of the brightness transform 
and is a measure of the spatial extent of the source relative to the resolution 
of the Interferometer for the given observation. It Is equal to the ratio of 
"correlated flux" to total flux and is defined by 

Vj(u,v,(d) h 


where S_((o) is the total power emitted by all parts of the source at frequency u. 
For a point source, the fringe visibility is unity for all baselines. For a dif- 
fuse source, the visibility is unity for a "null" baseline and decreases on 
average as baseline length increases. The decrease can become substantial for 
baselines satisfying 


RjujV,(i»2 

R(0,0,a>) 


|R(u,v,w)| 

Sp(a.) 


(A30) 

(A31) 




(A32) 


where dg is the source diameter in radians and X is the radio wavelength. 

In addition to this normalization of the brightness transform, we must 
normalize the voltage signal in Eq. (A4). Using Bqs. (All) and A14), one can 
easily show that expectation of the squared voltage is given by 



f fo 


) du dfl 


(A33) 


(A34) 


where we have again neglected the weak dependence of y< on k and evaluate y^ 
at the reference position Using Eqs. (23), (A6), iA22), (A31) and (A34), we 
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can now conpute Che croas-correlaclon fwictlon for Che normalized volcage conpo- 
nencs for Che nacural source: 



(AJ5) 


(A36) 


where 


5 / Sp(u)) G^(?^)du) / Sp(o)) G^(y^)du, 


(A37) 


and, like Eq. 27, 


(b- ■ U) (t* + T + T + T - T + T - lb) 
f g b I a ml mj 


* *h - "hi ’ml * *8 


(A38) 


in which Tb> ti. Tat and are differences between sCaCion quancicies In the 
order J minus i. 

As explained in the text, the quantity ~ difference between 

the errors in the bit times of the two stations, we have pulled out of the 
Integral, recognizing that resolution of the source will vary very little across 
a 2-HHz passband located at RF. A phase term has been added to fringe phase 
to represent the phase term from the brightness transform R. Further, we have 
used the definition of the instrumental phase from Eq. (A6), properly offset in 
time. In analogy with Eq. (26), a definition of an interim geometric delay has 
been made: 


o)T ’ H k • B 
gar 


• k. • - *•)) 


(A39) 

(AAO) 


where we have used the reference position k^ of the ext«ided source. 
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As iadlcsted in Bqs. (16) and (18), cIm final expression (Bq. 19) for the 
cross-correlation function of an extended source is obtained b> ':"Uiplping 
Bq. A36 by 2/« rad by a noise ta^ierature factor. One can easily snow tlMt this 
result for an extended source is the sane as that for a point source except that 
(a) an nplitude fact.;);' (vf 1 1> aust be included, (b) a phase shift due to 
structure auat be added to the phasq, and (c) the intaria delay is coaputed 
relative to the reference position ka« 
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APPENDIX B 

BIT MANIPULATIONS IN THE CIT/JPL MARK II CORRELATOR 


The JPL/CIT Mark II correlator has the capability of simultaneously 
processing up to five stations, including all 10 possible baselines. There are 
two stages in processing: single-station delay offsetting followed by two-station 

(pairwise) operations. Figure B1 illustrates these two types of operations in a 
5-statlon application and shows the numbering system for station pairs. Fig- 
ure B2 presents in greater detail the operations performed on a particular pair. 
Each bitstream is offset by a single-station model delay whose geometric compo- 
nent is the transit time between the station and the earth center of mass. The 
single-station delay for station 1 is given by 




ml 


1..4 + T 4 

ti cl 


(Bl) 


where the first term is the geometric delay relative to the earth center of « 
mass, expressed in terms of station location Xgii and apparent wave direction kn 
(both in geocentric coordinates). The term Tti represents the troposphere delay 
along the raypath to antenna 1, while represents the station clock error 
(offset and rate). The single-station approach works because the resulting dif- 
ferential offset between the two members of a pair of bitstreams will contain the 
differences for all important delay components. For example, the resulting dif- 
ference in geometric delays can easily be shown to equal the desired geometric 
delay (Eq. 26) between antennas (except for small effects due to argument offsets 
that are accounted for in post-correlatiou software). 

The pairwise operations illustrated in Fig. B2 consist of the following 
steps. Since the difference between two rounded numbers is not necessarily equal 
to the rounded difference between the raw numbers, the relative offset between 
two separately rounded and offset bitstreams can be in error by as much as ±1 bit 
relative to the desired rounded two-station delay. Thus a pairwise vernier cor- 
rection (1, 2, or 3 bits) is applied to the first bitstream of each station pair, 
where a 2-blt delay corresponds to zero correction. The same bitstream is then 
simultaneously subjected to each of a number (16) of lags. At the same time, the 
second bitstream is delayed by 10 bits so that its pairwise delay will match the 
total pairwise delay (for zero vernier correction) applied to the first bitstream 
in order to reach the reference lag (#8) for that bitstream (i.e., 2-blt vernier 
+ 8 lags = 10 bits). The second bitstream is then multiplied in quadrature by two 
pairwise stopping functions (cos($>tg, sln(|>g,) with phase (and rate) given by 
Eq. 28. The computation of phase through Eq. (28) must be based on total delays 
(Tjni) used by the correlator in offsetting to the reference (central) lag. That 
is, for each bitstream, the total unquantlzed delay is equal to the single-station 
delay (Eq. Bl) plus 10 bits (i.e., x ■ Xnii + 10 bits). The quantized delay 
XQj^ in Eq. 30 is equal to rounded to the nearest bit except that for the 
first bitstream is shifted by (-1, 0, +1) bit according to the vernier delay 
setting. After these pairwise operations, the two bitstreams (one now complex) 
are multiplied bitwise and a complex sum over bits is performed for each of the 
16 lags applied to the first bitstream. 
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Fig. Bl. Block diagram of the CIT/JFL Mark II multistation/correlator 
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SINGLE-STATION OPERATIONS 



71 


fig. B2. Sch— tic diagrm of two-sutlon crow-correUtlon .1th the CII/JPL c«t.Utor 











Figure B3 illustrates in terms of devised bits these bit manipulations for ! 

two bitstreams as they pass through the correlator. For the two possible levels 
of the recorded signal, a plus sign (minus sign) is used to denote a -<-1 (-1). 

Contrary to the impression created by the figure, the bits are processed serially, 
not simultaneously. The figure depicts a correlation Interval of 8 bits even 
though correlation (averaging) Intervals are typically 0.2 - 1.0 seconds 
(0.8 - 4 megabits) in length. No signal pattern has been made evident in the 
two bitstreams since the signal is typically nnich less than system noise for each 
1-bit sample point. As indicated above, the first step is to offset both bit- 
streams by their respective single-station delays given by Eq. B1 but rounded to 
the nearest bit. For illustration, those delays have been set equal to very 
small numbers (<3 bits) even though in practice they could be as large as 
20 msec (8 x 10^ bits). After single-station offsetting, one bit-stream (left) ' 

is subjected to the 1-bit vernier correction for roundoff errors, while the other ; 

bitstream is multiplied by the quadrature stopping functions. For simplicity in | 
the pairwise operations, the vernier delay is represented as (-1, 0, +1) bit 
rather than (1, 2, 3) bit; only one lag offset (l.e., no offset) is shown for | 

the left bitstream; and, consequently, the related 10-bit delay that would have 
been required for referencing the right bitstream to the central lag (the* only 
lag In the example) Is unnecessary. The stopping sinusoids are quantized to ' 

three amplitude levels equal to (-1, 0, +1). Over a cycle of stopped phase, 
the quantization subintervals are 3/8 cycle at +1, 1/8 cycle at 0, 3/8 cycle at 
-1 and 1/8 cycle at 0, with appropriate shifting for sine and cosine. One can 
readily show that the amplitude of the first harmonic Is equal to 1.176 for this 
quantization scheme. For our example, the fringe rate has been set equal to 
(8 blts)~^ or 0.5 MHz in order to make each quantization subinterval equal to a 
small Integer number of bits. This rate, of course, is much greater than the 
normal VLBl rates, which are 25 kHz or less at X-band. After multiplication by 
the quadrature stopping functions, the second bitstream has branched into two 
"counter-rotated" bitstreams, with some bits blanked out by the zeros in the 
stopping function. 

The two bitstreams derived from the second station are separately multiplied ' 
by the first bitstream to give two cross-correlated bitstreams. The bits in each 
of the product bitstreams are summed over the correlation Interval to obtain the 
excess (or deficit) of +l's relative to -I's. When divided by the total number 
of bits in the sum (see in Eq. 30), the two sums become the real and imaginary 
parts of the complex stopped fringes for which the theoretical expectation value 
is given by Eq. 37 with a^ = 1.176 and £ = 0. The time tag (i.e., bit time 
tG = tj,i = tbj in Eq. 29) for the correlation coefficients is the center of the 
sum interval in (unshlfteay UT seconds. 

Phase calibrator tones are extracted in the correlator by means of an artl- | 

flcial application of the cross-correlation operations. To set up for single- j 

station tone stopping, each playback bitstream is cross-correlated with (multiplied i 
by) an artificial bitstream filled with +l’s. The artificial "DC" bitstreams | 

are constructed in some or all of the bitstream channels left unoccupied by real j 

data. The model phase used for quadrature tone stopping is made to equal Eq. 60 | 

by an artificial assignment of the parameter values in the interferometer phase 
model in Eq. (28). Since tone stopping requires the use of unoccupied bitstream ) 
channels (see Fig. Bl), the maximum number of tones/frequency channel that can be | 
^ extracted in a single pass through the correlator (without time-multiplexing 

I between tones within a channel) is 5 - N^ is the number of real bitstreams being 

i 

I 
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Fig. B3. Bit manipulations with the CIT/JPL Mark II correlator 


processed. Thus, for two-antenna processing, a naximum of three continuous 
tones/frequency-channel can be obtained. By correlator design, more tones/ 
channel could be obtained by time-multiplexing between tone phase models. How- 
ever, as one increases the number of input bitstreams, increases the number of 
time-multiplexed frequency channels, decreases the dwell time per channel or 
increases the number of time-multiplexed tones within a channel, greater demand 
is placed on the available monory and computing power of the correlator. It is 
beyond the scope of this report to discuss correlator limits for these 
parameters. 
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APPENDIX C 

REFORMULATION AND ILLUSTRATION OF FRINGES IN THE LAG DOMAIN 


This appendix reforaulates the stopped fringes in Section VI to place them 
in the form that is commonly used in the lag domain. An example of stopped 
fringes is presented. 


The identity 


iioAi . -wAt , , 

e i e £ = 1 


(Cl) 


can be used to place the stopped fringes of Eq. (37) in the form 


<u,(tg» ■ .J -f e 


llAi|ij-((i)-<Dj^^)£b] 


(C2) 


where the delay function is given by 




s i j S_(») G^(y^) Cj(,j 


) 


l(ta)-U))AT, 


du 


(C3) 


and where the overall residual fringe phase at the reference lag (£-0) is given 
by 


= A(|»j^ + u)At + ♦j + •J'g 


(C4) 


and residual delay by 


At^ » At - £b (C5) 

The frequency u> will be equ^l to the centroid of the bandpass product. One can 
show that this choice for u will relegate virtually all of the time dependence 
in the fringes to the fringe phase \l>f. 


Suppose that the natural source spectrum Sp(u>) is flat and that the bandpass 
shapes at the two stations are Identical. Under these assumptions, the delay 
function becomes 


D (At) 
c 




/ 


G(yj)G(yj)F^e^^“"“^^‘' do) 


(C6) 
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where In Eq. (21) becomes 


'N 


/ 


G ((A))du) 


(C7) 


Given the bandpass shape, the quantization filter Fq, and the doppler shift in 
y£ and y j , this Integral can be readily evaluated. The sinq>le case of a 
rectangular bandpass of width W{) proceeds as follows. The BLKO system operates 
in the single-sideband mode in which the total mixing frequency is placed at the 
edge of the 2-MHz passband. If the doppler shift at station J is Dj and the 
total mixing frequency is f j , then the centroid frequency will be 


0) 

2ir 


2 


(C8) 


where the plus (minus) sign is used for upper (lower) sideband. Effective 
bandpass width will restrict the limits on the integral in Eq. (C7) to the 
range 


0) = u) ± Wj/2 

d 


(C9) 


where W^j is the bandpass overlap after doppler-shif ting given by 

provided there is overlap. Given these values, the integral in Eq. C6 becomes 


D (At) 
c 


W. 8in(iTWjAT) 
d d 

Wb itW^At 


(Cll) 


where we have assumed Fq = 1. The factor which lies between zero and one, 

is the ratio of actual overlap to total bandwidth. It gives the fractional 
power remaining after doppler shifting has caused misalignment of the two station 
bandpasses and has thereby dropped power at the bandpass edges. For the BLKO 
system, such losses are almost negligible since doppler shifting is very small 
relative to a 2-MHz bandwidth. For example, the maximum doppler shift at X-band 
is about 20 kHz so that Wj/Wb 1: 0.99. (In the BLKO system, f^ - fj is often zero 
since the same mixing frequencies are often used at both stations. In narrowband 
systems, the mixing frequencies are sometimes given different values so that their 
difference compensates for the bandpass misalignment caused by doppler shifting.) 


Thus for a rectangular bandpass, the stopped fringes are given by 


r W. 8in(irW.AT.) i(Ai(»,-(u)-a). . )ib] 
m d d % X hi 

" ®1 IT * 


(C12) 
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where Ai|/£ Is given by Eq. (C4). If power loss due to doppler shifting is 
negligible (W^/Wf, ■ 1), the maximum amplitude of the fringes is given by 
airg|/ir for perfect bitstream alignment (At * 0). 

Examples of lag-domain fringes that correspond to the frequency-domain 
fringes in Fig. 3 are given in Fig. Cl and are plotted as a function of time 
for three lags (i). Note the phase shift of about ir/2 between the sinusoids for 
adjacent lags. This shift is a consequence of the term (u-ti>||j^)ib in Bq. (C2), 
which is caused by lag £ffsettlng of bitstream i. Since the average (centroid) 
frequency at baseband (o> - is about 1 MHz, a shift of one lag will produce 

a phase shift equal to 


2^ (“-Whi)b » ±1 MHz * 1 bit (C13) 

* ±0.25 cycle 

where the ± pertains to upper and lower sidebands. The actual phase shift is 
slightly smaller (*0.22 cycles) since the bandpass falls off before reaching 
2 MHz. In fringe stopping, the same stopping phase is used for the fringes 
for all lag offsets so that this phase jump between lags is not removed. Also 
displayed in this figure are fit fringes based on the model in Eq. (C12), where 
the bandwidth (Wjj=Wj) has been set equal to 1.8 MHz. In the fit, the behavior 
of the amplitude as a function of lag allows determination of the fringe ampli- 
tude and BSA delay. For the first section, one can visually fit a sin x/x func- 
tion to the amplitude vs lag in Fig. Cl and verify the fringe amplitude of -800 
and the BSA delay of -0.3 bit that result from phase-tracking. By noting the zero 
crossings, one can also verify the fringe rate of about 21 mHz and the fringe 
phase of about 0.5 cycle at the indicated midsection reference time. The SNR 
per fringe point is about 8 in this example. 

With regard to amplitude, the solve-for amplitude in the lag-domain should 
be approximately equal to a^ (see Eq. C2) since no amplitude scaling is 

provided in the lag-domain model fringes. This amplitude is smaller than the 
frequency-domain result by about a factor of 2.0 due to amplitude scaling through 
BAMPL, as explained in Section X. 
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SECTION I 


SECTION 2 


A^.809 ±24, 4,-21.30 t0.27mHs 

. 0.497 t0.0048 cycN, r, - *0.321 10.054 lap 


A^.633 ±24, ^,-24.00 ±0.27 «Hi 

4, • 1 .890 ±0.0046 eyeb, r, - -0.321 ±0.053 lap 



TIME, wc 


Fig, Cl. Long-baseline example of stopped fringes in the lag domain 
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APPENDIX D 


QUANTIZATION FILTER 


This appendix obtains an expression for the quantization (fractlonal-blt- 
shlft) filter by evaluating the sum 


Fq(z) • (Dl) 

‘ 'g 

Ue will assume in the following analysis that the corr<>lator keeps the quantized 
model delay within 0.5 bit of the model delay (i.e., |Tg, - Tq| 1 b/2). If one 
assumes that delay rate (t ) is constant across the sum (correlation) interval 
(which It very nearly will oe for correlation intervals as short as 1 sec or 
less), then the sum over time can be converted to an integral over t. The 
integral over t, however, must be broken into intervals with boundaries at bit 
Jump points. 


We will model the general case in which the quantization error 
(At ■ Tg, - Tg,) exhibits the following behavior. In the first interval with 
continuous behavior, let At progress from At]^ to b/2. This will be followed 
by Ng intervals in which At traverses its full range from -b/2 to b/2. In the 
final interval. At then changes from -b/2 to At 2 . By evaluating the appropriate 
integral in each interval and summing, one can show that the quantization filter 
will be equal to 


F,(I) 



IzAt^ IzATj 
e -e 


Iz 


+ Njb 


sin zb/2 
zb/2 


(D2) 


where t^^^^ and are the initial and final values of Tgj,. respectively, and 

where Nj is the number of "Jumps" in the quantized delay Tg, over the full 

correlation Interval (Nj - N^ + 1 ■ Tg ,2 - ^mi)* Equation D2 is quite general 
and is valid when there are no Jumps in the quantized delay (Nj ■ 0) and when 
the delay rate is negative (Tg,2 “ '^ml ^ ® 

Two limiting cases are of Interest. First, in the limit of large Nj, the 
quantization filter introduces no phase shifts and, as one would expect, becomes 
a pure sin x/x filter. This is also true whenever Atj ■ At 2 ■ 0 and |Nj| > 0. 

Second, as the delay rate approaches zero, we will have Tg ,2 - Tg,]^ -*• 0 and 

Kj ■ 0. In this limit, we obtain the expected result 

F (z) - exp(lz(T - T )] (D3) 

q mm 

where t Is the delay at the middle of the correlation interval, 
m 
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The most important feature of the foniulation of the quantisation filter in 
Eq. (D2) is that only the initial and final values for Tq are needed to compute 
the composite effect of quantization over the whole interval. Given the initial 
and final values, one can easily coiq>ute Tq 2 ” '^ml* Nj. For this 

reason, Tmi and 1^2 values for each correlation Interval are numerically passed 
to the modeling subroutine In PMASOR. 
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APPENDIX E 


ERRORS IN THE LAG-TO-FRBQUENCY TRANSFORMATION 


In the transfona&tlon of the fringes from the lag domain to the frequency 
domain, the use of a limited number of lags can introduce errors in the output 
phase and amplitude. This appendix develops a model to assess the magnitude 
of such errors as a function of the total number of lags used in the transform. 

To stress the important features in the lag transform, we cast the fringes 
of Eq. (37) in the form 


<u^> « 


H(a>')e 


ioi'(AT'-ib) 


where H(u') is equal to the bandpass product SpC^G, and oi* is the translated 
frequency u* > We will neglect doppler shitting and will assume that the 

two bandpasses are identical (G^Gj “ G^). The delay At* is the residual after all 
model terms are subtracted. In this appendix, we will let the integer £ equal 
zero at the central lag, which is chosen, of course, to cancel all but the frac- 
tional bit part of At*. Thus the "remaining" At* will fall between ±0.5 bit. 
Although we will not do it here, it can be shown that the quantization filter Fq 
has little effect on the transformation errors and can be neglected. All phase 
terms independent of frequency have been discarded since they will only cause a 
phase shift in the transformed fringes and have no effect on lag transform errors 
(although they are otherwise important). 

The error calculation proceeds as follows. For a given bandpass-product 
shape and delay At', the Integral in Eq. (El) is calculated numerically for each 
of N^ values of £ centered about the central lag, thereby simulating the proces- 
sor output. The resulting N|^ complex numbers are then subjected to the limited 
inverse transformation in Fq. (40), which is the way the postcorrelation software 
transforms the fringes to the frequency domain. The resulting complex numbers 
in the frequency domain will belong to frequency "bins" centered at the frequen- 
cies in Eq. (43). These frequency-domain values should be approximately equal to 
the complex input Integrand in Eq. (El) (without £s) at each frequency bin. 

An example of the process is shown in Figs. Ela and Elb for a nine-lag trans- 
form and a Butterworth filter approximating the BLKO bandpass. The solid lines 
represent the input amplitude (G^) and phase ((d'At') in the frequency domain, 
while the points represent the limited-lag transform output that shfmld be 
approximately equal to this input. For this particular case, the output aiiq>ll- 
tude deviates from the input amplitude by as much as 9X, relative to ma x im u m ampli- 
tude, while the output phase oscillates about the input phase by about ±0.5 milll- 
cycles. The directions of the deviations in this example are fairly typical, 
although the m^nitudes are not. 
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Pig. Ela. An exaiq;>le of lintted-lag transfornatlon errors 
for an ideal Butcerworth filter 
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Fig. Elb. Plot of phase errors in Fig. Ela 
on a finer scale 
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Since complete results for all cases would be too lengthy, we have 
suDunarlzed in Table El the transformation errors calculated for certain useful 
cases. Bandpass A roughly approximates the iu.kO bandpass, while bandpass B is 
somewhat wider as explained in the table caption. The table was constructed 
as follows. For a given bandpass shape and given lag count, the transformation 
errors were calculated for various values of At' between ±0.5 bit. Results are 
shown only for odd bln frequencies [n/(2N£), n odd]. It was found that the 
amplitude errors for fixed bandpass and fixed lag count (Nj^) are virtually inde- 
pendent of At ' while the phase errors depend strongly on At ' . The table gives 
only the maximum error in bin phase across the bandpass for the worst case of 
At*. In addition, the table lists the worst-case error in aggregate phase 
(referenced to the bandpass center) obtained from a weighted least-squares fit 
of a linear function of frequency to the bln phase errors. In the fit, the bin- 
phase values were weighted by the square of the bln amplitude. (As suggested by 
Section X, such a fit is implicitly performed when the transformed fringes are 
phase- tracked. ) Due to the oscillation of the bin errors, the worst-case error 
in aggregate phase is typically much smaller than the bin phase errors. The error 
in aggregate phase falls to an acceptable level (^0.001 cycle) at 7 or 8 lags. 

From the table, it appears that a 9-lag transform yields the smallest phase 
errors, both bin and aggregate, for Nj^ below 12. For » 9, bin-phase errors 
are less than a mlllicycle and aggregate phase error is less than 0.1 millicycle. 

With regard to amplitude errors, the table also gives, for each case, (a) 
the largest discrepancy between input amplitude and output amplitude across the 
passband, expressed as a percentage of input amplitude, and (b) an estimate of 
aggregate amplitude error for the case in which bln amplitude in fringe fitting 
is approximated with the actual bandpass amplitude in Eq. (71) (l.e. , the bin 
amplitude for -> <»). From the table, we see that the worst bin-amplitude error 
usually occurs at the first or last bln in the bandpass, where edge effects are 
most significant. The aggregate amplitude error is 2 or 3% when the lag number 
is 9. For very accurate amplitude work, this error can be eliminated in principle 
by making accurate measurements of the bandpass shape and calculating bln ampli- 
tudes by means of the simulation analysis described above. Although it does not 
appear to be necessary for any current applications, bin-phase errors could be 
largely eliminated in a similar fashion. 

A word of caution is necessary in using these results. While the numerical 
values can provide a useful rough estimate of transformation errors, the simulated 
errors might deviate considerably from actual errors. The main reason is that 
actual bandpass shapes can deviate considerably from the assumed ideal Butterworth 
shapes. More accurate representations of bandpass shape might lead to large 
changes in aggregate phase errors, which were reduced to small values by cancel- 
lation of large bin-phase errors of varying sign. 

The shape of the transform bandpass obtained in this appendix can be used 
to compute nominal values for BAMPL, the four-element array in PHASOR for nor- 
malizing bln amplitudes in frequency-domain fits. On the basis of Eqs. (45) 
and (49), one can readily show that 

BAMPL = 1-176 (0.94, 1.0, 1.02, 0.89) (E2) 
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Table El. Transformation errors from limited- lag transforms 
[Bandpass A (sBLKD) is an ideal 13-pole 
Butterworth filter with -3 dB at 1.76 MHz; 
Bandpass B is an ideal 7-pole Butterworth 
with -3 dB at 1.83 MHz] 


Lag 

count 

Ni 

Band 

pass 

Largest 

bin-phase 

error, 

mcycle 

Largest 

composite 

phase 

error, 

mcycle 

Largest 

bin 

amplitude 

error,* 

% 

Composite 

amplitude 

error* 

% 

4 

A 

10.4 

1.1 

-21 

-12 


B 

6.0 

0.5 

-13 

-8 

5 

A 

5.8 

2.5 

-7 

-2.0 


B 

5.7 

2.4 

-7 

-2.1 

6 

A 

10.3 

1.0 

-17 

-5.8 


B 

5.1 

0.8 

-9 

-3.7 

7 

A 

1.6 

0.5 

-7 

-3.3 


B 

2.0 

0.6 

-7 

-2.6 

8 

A 

9.0 

0.25 

+5 

-0.5 


B 

3.7 

0.03 

-6 

-1.5 

9 

A 

0.55 

0.05 

-9 

-3.1 


B 

0.51 

0.05 

-6 

-1.8 

10 

A 

7.5 

0.02 

-7 

-0.2 


B 

2.4 

0.10 

-7 

-1.2 

11 

A 

1.7 

0.17 

-10 

-2.5 


B 

0.2 

0.02 

-7 

-1.4 

12 

A 

5.3 

0.15 

+7 

-0.7 


B 

1.3 

0.02 

-6 

-1.0 


*As percent of maximum input amplitude. 
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where the first factor Is ax and 0.88 is the value of Ga in Eq. (50) computed on 
the basis of the bandpass shape (solid curve) in Fig. Ela. (More exactly, for 
constant Sp and identical band passes at both stations, Ga becomes after Sp 
has cancelled in niraerator and denominator of Eq. 49). The four numbers in tne 
array are the four bin amplitudes in Fig. Ela. The new factor 0.75 enters because 
the correlator divides by the number of unblanked bits (0.75N^) rather than the 
total number of bits (N^) . The overall value for the BAMPL array becomes 

BAMPL = (1.07, 1.14, 1.16, 1.01) (E3) 

These values for BAMPL are used as the default values in PHASOR. These BAMPL 
values correspond to a "b-factor" of about 2.0 whereas experimental determinations 
of the b-factor have resulted in a value about 25Z larger. That is, if the 
default BAMPL array is used, the output amplitudes are reported to be about 20% 
too small. The cause of the bulk of this discrepancy is not currently known. 


>r. 
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APPENDIX F 


THE PHASE CALIBRATION SIGNAL 


This appendix treats the phase calibrator signal and Includes an analysis 
of (a) tone power, (b) the effect of bilevel clipping on the recorded signal, 
(c) tone stopping, and (d) stopped-tone SNR. Since It Is not necessary for 
these topics, this appendix will not decompose tone phase Into signal-path 
components but will treat It as a single variable. An analysis of the compo- 
nents of tone phase Is given In the text. 

The analog baseband signal (before clipping) can be represented In the 

form 


V(t) = S(t) + Ojjd(t) (FI) 

where S(t) Is the total calibrator signal and the second term Is additive noise, 
which Includes Instrumental effects, sky noise and source noise. We have 
normalized the noise to unity (l.e>, <d^> » 1) and Included a factor equal 
to the rms voltage of the noise term. 

The calibrator signal can be decomposed In terms of Its harmonics to give 


V(t) 



+ 0jj8(t) 


(F2) 


where Is the phase of the n^^ tone and Vj, Is the maximum voltage of the n^^ 
tone. For more detail concerning 41 ,^, see Eq. (57) and the discussion In the 
text. The amplitude v.^ will be discussed below. 

Total power of the calibrator signal Is given by 


P 

s 



while the power of the noise term Is given by 

„ 2 

P <* o 
N N 


(F3) 


(F4) 


If we specify total calibrator power relative to total noise power (both within 
passband), we can express the tone amplitude In terms of 0 }|. That Is, If the 
ratio (Ps/^N^ total calibrator power to total noise Is e, then the ratio (SNR) 
of tone amplitude to rms noise voltage becomes 


V 





(F5) 
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where is the number of Cones In Che passband ,_and v is maximum tone 
ampliCuae, ideally assumed to be the same (vn » v) for all tones in the passband. 
(It is unnecessary for our purposes to resolve the uncertain status of tones 
on the passband edges.) We currently plan to keep the total calibrator power 
low — about 2 % of the system noise. Thus, if there are three tones in the pass- 
band, the SNR for a tone becomes 


V 


I 


2 X 0.02 


0.12 


(F6) 


(For all examples, we will assume a nominal system with a 2 % power level and 
3 Cones/ passband. It is outside the scope of this report Co determine optimum 
values for these parameters.) This "single-sample" SNR is large compared to 
that of typical Intercontinental cross-correlation fringes, which are usually 
in the range of 0.001 to 0.04. The tone SNR will be modified slightly below 
when two-level sampling is taken into account and will Increase when many bits are 
collected in the tone-stopping process. 


The maximum amplitude of the signal S occurs whenever all tones are in 
phase and is given by 


S 

max 


N, V 
n 


(F7) 


This maximum occurs, of course, every time a calibrator pulse reaches the 
recorder and equals about 0.35 for the example above. (The pulse width at 
baseband is approximately equal to the inverse of Che system bandwidth when a 
number of tones are in the passband.) 


As indicated in the text, the baseband signal is subjected Co bilevel 
sampling and recorded at 4 Mbit/sec. We will use a tilde (~) to denote Che 
bilevel value (±1) at each sample point. Unlike interferometric cross- 
correlation in which both the signal and noise are random, here the tone signal 
is deterministic and the noise d(t) is random. Thus the expectation of a 
sample point becomes 


<V> 



Q(V)P(V)dV 


where the sampling function is given by 

Q(V) = +1, V > 0 


= -1, V < 0 


(F8) 


(F9) 
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and where the probability distribution of V (about S in Eq. FI) is given by 


P(V) 


1 expt-(V-S)^/(20jj^)J 

ir o„ 


(FIO) 


Equation (F 8 ) becomes 


so that 



(Fll) 


(F12) 


(F13) 


if one neglects higher order terms. As indicated above, the maximum value of 
S/ 0 {{ in our example will be about 0.3S so that the maximum of the cubic term will 
be about 2 % relative to the first term (before tone stopping). 

We are now prepared to analyze tone-stopping in the correlator. Let V]^ be 
the recorded signal at time tj^. For a given station, the correlator will 
counter-rotate (tone-stop) a given tone as follows 


''.n - \ 

' k=l 

where is the model phase for the n*-^ tone at time t)^ and is the number 
of sample points In the correlation Interval. We will not take into account the 
fact that the BLKO correlator uses a trllevel quantized model for the stopping 
sinusoids. Such trllevel sinusoids possess higher-order harmonics that can cause 
problems if tones are not properly placed in the passband. A discussion of these 
problems is outside the SC 07 .! of this report. Trllevel quantization also 
results in a slight decrease In SNR and a change in stopped-tone amplitude, but 
these effects are not significant in most applications. However, if absolute 
calibration of tone amplitude Is desired, the amplitude effect must be 
considered. 


89 


Given Eq. (F13), the expectation value of the stopped tone given in 
Eq. (F14) can be calculated as follows: 


N. 


<V > 
sn 


k-l 


<V^> 


^ k=l “ 


(F15) 


(F16) 


where S|^ is the total tone signal S at time t|(. The neglected cubic term in 
Eq. (F12) can be represented as a sum of three-tone beat notes. One can show 
that during tone stopping a given three-tone beat note will be reduced to an 
insignificant level provided it possesses a different frequency than the original 
tone frequencies. If this £ondltlon is not satisfied, the offending three-tone 
beat note in a system with v/ou = 0.12 could fractionally contribute in ampli- 
tude as much as 


1 

4 


* 



0.0006 


(F17) 


or 0.06% to the tone under consideration. Although this is not much by itself, 
many equally spaced tones can generate many three-tone beat notes near a ^iven 
original tone frequency and the combined effect might be non-negliglble. This 
problem will deserve more thorough study if many tones are ever used. 

Based on the decomposition of S in the first term of Eq. (F2), the stopped 
tone in Eq. (F16) becomes 



(F18) 


if we assume that the sum note and the other "high frequency" tones sum to 
negligible levels. The phase difference 4>n - 'I'n* presumed to be nearly constant, 
is the phase difference at the middle of the sum interval. (Ue will assume stopped- 
tone frequency is negligibly small. If it is not, there are straightforward, 
simple ways to overcome the difficulty, such as the procedure used in the PCAL 
program described in Appendix G.) 
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(F19) 


The syston noise on the stopped tone can be calculated as follows: 

V y sn sn / *^ 

P 

where is the real part of the stopped tone in Eq. ('^lA). Using Eq. (F14) one 
obtains 


o 


2 

V 


■ “''k" "''‘j 

t kt *• 


These expectation values can be evaluated to obtain 


(F20) 


”5 ' 2T * f H “<’> ‘=‘” 2* 

TjtO 

where the sum is over all bits except t ■ 0, R(t) is the bitstream autocorrela- 
tion function for the analog signals and Vn is the baseband frequency of the n^^ 
tone. This derivation assumes that Interblt correlations of system noise are 
small. The first term represents the noise that would be present If there were 
no correlation between bits. The second term accounts for small correlations 
between bits and will be about 10% or less, the exact size depending on the 
shape of the bandpass, the sample rate relative to bandwidth, and the value of 
the baseband tone frequency Vn. Note that, if the small second term can be 
neglected, the noise on a stopped tone Is the same for all tones. 

Stopped-tone SNR will be defined as the maximum amplitude (modulus of 
Eq. 18) of the stopped tone divided by the rms noise: 


SNR 


1 _v_ X 

JTiT 


which through the use of Eqs. F5 and F21 becomes 





bWK • ^ 

wNh 


(F22) 


(F23) 


where we have neglected the small Interbit-correlation term in Eq. F21. In analogy 
with fringe phase calculations, the system noise error In the stopped-tone will be 

a, - SNR"^ (SNR » 1) (F24) 

•P 


91 



The following example will help Illustrate the size of the system noise 
error in phase calibrator corrections. Again, suppose we have a 2Z total power 
level for the calibrator signal with 3 tones across the passband. When 4 seconds 
of BLKO data (record rate > 4 Mbit /sec) are reduced for a given tone, we obtain 
from EQ. (F23) 


SNR 


1 2 X 0.02 X 4 

■yj 3it 


X 4 X 10 


260 


(F25) 


The associated phase error will be 


% * Ym * 0.0006 cycle 


(F26) 


This phase error Is the system noise error in the phase from one tone. 
When fringe phase is calibrated, a number of tones (Ny) in the passband will 
probably be used, in which case the system noise error in the overall phase 
would be 


T - 

0 . » 


JL 


JX 


jr h 

2N,e N 
t u 


(F27) 


For our example, the overall phase error becomes 

at 0.00035 cycle 


(F28) 


If all 3 tones are used. For a 40-MHz spanned bandwidth in a bandwidth s 3 mthesis 
delay measurement, the corresponding delay error will be 


^ /2 X 0.00035 cycle 
~ 40 MHz 


0.012 nsec 


0.36 cm 


(F29) 


The factor of /2 is a result of the combination of two channels in the BWS 
process. 

Before closing the discussion of the calibrator signal, a more detailed 
description of tone amplitude is In order. If the ideal calibrator signal in 
Fig. 2 is decomposed into Its harmonics, it can be shown that the amplitude of 
the n^h harmonic or tone will be given by 


sin((u T /2) 
F (o) ) - C 

P " 


(F30) 
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where tone frequency utn Is a multiple (nup) of the nominal pulse rate u. and 
C is a constant determined by input power. With a pulse width of Tp « 20 psec, 
the tone amplitude at X-band {8.3 GHs) will fall to about 

(u) T /2)^ 

F (w ) « 1 ^ - 0.95 (F31) 

p X 0 

which represents a loss of about 5% relative to the maximum amplitude at lower 
frequencies. The amplitude change across the passband at S- or X-band will be 
even smaller. For example, if the passband is 100 MHz, the amplitude change 
will be about 0.1% across the band at X-band. Thus, changes in tone amplitude 
across the passband due to pulse shape appear to be negligible for an ideal 
tone generator. Actual tone amplitudes might deviate considerably from the 
ideal behavior described here. 
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APPENDIX G 

THE CURRENT METHOD FOR PHASE CALIBRATION 


In Che text. It Is assumed that phase calibration is applied at an early 
stage in postcorrelation processing and that the correction is applied to 
t'ue stopped fringes in the frequency domain by means of a frequency-specific 
counter-rotation of the phase in each correlation interval. However, the cur- 
rent BLKO system postpones phase calibration until a later st^e of post- 
correlation processing at which point a program called MERGE applies corrections 
to the composite values for phase, BWS delay, BSA delay and phase-delay rate 
that are passed by PHASOR. This version of phase calibration was coded first 
because it could be more easily incorporated into existing software. The compo- 
site interferometer observables are extracted as described in Section X except that 
the phase-tracked fringes are not subjected to phase calibration. A block diagram 
of the current data reduction steps is shown in Fig. Gl. 

To obtain tone phase, the stopped tones produced by the correlator are 
passed to a special phase- tracking program (PCAL). In close analogy with PHASOR, 
this program fits the stopped-tone signal in separate time intervals (sections) 
within an observation with a model in which the three solve-for parameters are 
amplitude, phase and phase rate. As in PHASOR, the solve-for parameters 
obtained for the various sections are combined to obtain for each parameter a 
composite value for the whole observation. For each tone, the composite ampli- 
tude, phase and phase rate are specified at a reference time and frequency near 
the center of the observation. 

In general, tone phase cannot be subtracted from fringe phase directly 
because the reference times and frequencies of those phase values are different. 
First, tone phase is projected to the correct time using its measured phase rate; 
then calibration phase at the correct frequency is obtained by interpolating 
between the (time-corrected) phases of two tones. Typically, there are tones 
within ±0.5 MHz of the bandpass center. The resulting calibration phases are then 
subtracted from the composite fringe phase channel by channel to obtain a value 
for corrected fringe phase given theoretically by Eq. (77). For each observable, 
it should be understood that "subtracting" the tone values actually means sub- 
tracting the appropriate difference between stations, as in Eq. (64). 

For phase-delay rate from a given channel, the composite tone phase-rate 
values are used to Interpolate to the reference* frequency, and the resulting 
value is subtracted from the composite interferometer phase-rate value. For BSA 
delay from a given channel, the composite tone phases from that channel must be 
analyzed to obtain the effective slope ("group delay") as a function of fre- 
quency. In the case of two tones per channel, the slope is simply computed as 
phase difference divided by frequency difference. As mentioned in Section VIII 
and as illustrated below, one must have initial crude values for the slope in 
each channel In order to resolve integer-cycle aiid>igulties in tone phase. The 
correction Is then carried out by subtracting this tone group delay from inter- 
ferometer BSA delay for that channel. Finally, for BWS delay, the time- and 
frequency-corrected calibration phase values from the channels are combined 
in pairs as in Eq. (81a) to form tone "BWS delays" which are then subtracted 
from corresponding interferometer BWS delays. 
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STATION I STATIONS STATIONS 



FINAL VALUES FOE TOTAL FHASEi 
TOTAL DELAY AMO TOTAL 

FHASE*OELAY HATE 


(a) MUTIBASaiNE AND MULTIFKEQUENCY 
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STOFFINO USING CIUOE OEUY MODEL 

M DETECT AND FNASC'TEACK 
X-COEttLATION SIGNAL; FOEM 
A M II GU OUS DNS DELAYS; IMFEOVE 
DCUY MODEL 

(c) FHIGE-TIACK TONES 

(d) FEEFOBM FHASE CAUMATIOH EEMOVE 
INTEGEE-CVCU AMHOWTIES SY 
SINGU-OESEEVATION EANOWIOTN 
SYNTHESIS 

(•) FEEUMINAIIY FIT TO EESIMAL 

OKttVAILES FOE SINOU-SASaiNb 
CHECK INTfiOa<CYCU AMEIOMTIE^ 
DATA QUAUTY ASSESSMENT 

(0 CHOOSE «KST' OESEEVAILB 
(CHANNEL, CHANNa>FAIE>| 
TMJFOSFHEEE AND CHM^ MKnCU 
CALWATIONS; EESTOEE MODEL OnAYi 
EDIT DATA 


Fig. Gl. Block diagram of present data reduction in 


:c BUU) system 
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An example of the instrumental phase measured with tones is given in 
Fig. G2 for the BtKO dystem. The figure plots the phase for 19 tones in a 
2-MRs bandpass recorded at DSS 43 for channel 1 in our example observation. 
Although many tones were processed in order to produce this figure, no more 
than three tones per passband can be processed for a two-antenna experiment If 
only a single pass is made through the correlator, as discussed in Appendix B. 

As indicated in the figure, a linear function of frequency has been subtracted 
fiom tone phase in order to reveal nonlinear effects introduced by the instru- 
mentation. The large nonlinear trend that remains is due primarily to the narrow- 
band IF- to- video converter. For this example, the total tone phase at the 
reference frequency (0.9 MHx) is equal to 1.285 cycle and is obtained by adding 
the Indicated linear function to the plotted value. The time extrapolation 
mentioned above was negligible for this particular example since the reference 
time for tone phase was very close to that for the composite interferometer 
phase. 


Another feature of phase calibration should be noted in this example. To 
reconstruct instrumental phase across the passband, one must resolve a slope 
ambiguity equal to n/Af, where n is an unknown integer and Af is the separation 
between the two closest extracted tones. For the present example, the ambiguity is 
10 usee so that only crude calibration of the instrumental delays discussed in Sec- 
tion VIII is necessary to remove the slope ambiguity. Calibrations by L. E. Young 
give approximate values^ for the three delays in Eq. (62): 


(a) the epoch difference in the recorder clock and the ^-cal clock “ 
T. 


bj 


- T . • 5 usee, 
cj 


(b) the delay ("cable delay") from the clock reference point (<^-cal clock) 
to the ^-cal injection point ■ tyj ml. 3 usec. 


(c) the delay from the injection point to the recorder ■ txj m 2.4 usec. 


These terms give a slope of “(t. . - t . -f t . + t_.) - -9 usec. We will 
ignore the correlator reference tine ts sl&ce, iH'^ this example, tg is an 
integer multiple of 1 second and all tone frequencies and the reference fre- 
quency are integer multiples of 1 hertz. Thus, the resulting products of 
time and frequency each contribute integer cycles at desired saiq>le points. 


^The offset between the recorder clock and the 4~cal clock, T), - Tc, is obtained 
by subtracting the offsets separately measured for those two clocks relative to 
the station 1 ppo. When is measured, one must make certain that the 4-caI 
clock is based on the particular subset of positive-going zero-crossings (at 
the reference point discussed in Section III) that are allowed to become pulses 
at the output of the tone generator. Since, in our example, every SOth zero- 
crossing of a 5-NHz signal is converted to a pulse, the separation between 
passed pulses is 10 usec. This pulse separation is rei.ated to the 10-usec slope 
ambiguity for tone phase, which le due to the ambiguity in time tags for the 
"pulse-generatliig" zero-crossings at the reference point. Once tiiose time tags 
have been assigned, can be determined and the slope ambiguity can be worked 
out, as will be illustrated. 
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which would Incrwnent the slope In steps of 10 )isec< As shown in Section IK* 
the t 3 term In the slope disappesrs in the difference between stations anyway. 

In addition to these nK)deled delays, the artificial method used in the correla** 
tor for tone extraction introduces another unmodeled delay of 2.5 usee, which 
brings the total slope to about **11 usee. This a priori estimate of slope 
clearly selects as the closest ambiguity in slope the value -11.221 ysec, given 
that the ambiguity spacing is 10 usee. (Although the unmodeled delay introduced 
by the correlator does impact single-station phase and slope determination, its 
effects cancel in the difference between stations and will therefore be ignored 
after this point.) 

Nonlinear trends as large as those in Fig. 62 can lead to significant errors 
in phase calibration. For the reference frequency approach presented above, the 
true composite phase effect obtained by integrating across the passband can differ 
considerably from the observed phase at the reference frequency, as suggested by 
the asymmetrical weighting imposed on the large nonlinear phase excursions at the 
bandpass edges. Further, the slope (BSA delay correction) obtained by using two 
points in the central part of the bandpass will not be as steep as the actual 
"average” slope of phase. However, the actual errors in composite phase and 
slope are not as bad as this example suggests since an Interferometer forms the 
difference in phase between antennas and the above single antenna errors will 
cancel to the extent that the system bandpasses of the stations are identical. 

For BWS delay, further cancellation can take place between channels at the same 
station. It is beyond the scope of this report to assess the magnitude of 
phase calibration errors. The next generation system, the BLOCK II system, will 
reduce such errors through the use of more tones. 

Although the specific calculations will not be presented here, the composite 
tone phase extracted for channel 3 at DSS 43 for the same observation was equal 
to 1.007 cycle at the reference frequency. The results for the two channels can 
now be combined to obtain the DSS 43 tone correction for BWS delay for the (1,3) 
channel pair: 


At 


(1.285 - 1.007) cycle 
(2276 - 2271) MHz 


55.6 nsec 


(Gl) 


The same calculation can be carried out for the other station (DSS 14 in the 
example observation) to yield a BWS delay correction for that antenna. The 
results for the two antennas will be used in the next section in an Illustration 
of single-observation bandwidth synthesis. 
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PISCEDINQ 

APPffiDIX H 

81N6U>0B8BRVAI10N BANDHUira SYRraSSlS 


One of the goals of VLB! measurements is the precise determination of 
delay (i.«. • dTf in Eq. 77). One can easily show that syston noise error in 
delay is inversely proportional to the highest frequmcy (for phase-^elay) nr 
frequoacy spread (for BUS delay) in the observed spectrum* provided int^er- 
cycle ambiguities in the associated phase or phase difference can be tauoved. 

The BSA delay Tg is based on the continuous single-channel bandpass and becomes 
an unambiguous estimate of At£. However, since this bandpass is relatively 
narrow (2 MHz), the BSA delay is relatively Imprecise. By conqiarison, the phase 
delay (^s/'^z)* which is at RF, is an eztrmnely precise measurmnent of delay except 
for an integer-cycle andiigulty. One method for covering the range between these 
two extremes in frequency without integer-cycle ambiguities is bandwidth synthe- 
sis. In the BLRO bandwidth synthesis technique, one measures interferometer 
phase at selected frequencies (l.e., channels at wg “ Wg, at], •••) and then 
iteratively rmnoves the ambiguities in the difference phase between successively 
larger channel separations. 

This appendix describes single-observation bandwidth synthesis, rather than 
multi-observation synthesis as outlined by Purcell (Ref. 6). When phase calibra- 
tion is available, as assumed In this report, one can resolve the antblguities in 
phase by using only the information found within each observation. Without phase 
calibration, one must use a multi-observation approach that removes ambiguities 
only in a relative sense by comparing adjacent observations. The single-observation 
approach with phase calibration is preferable since it is sioqtler and faster and 
produces an absolute measurement of delay that allows clock synchronization 
measurements. 

The orig inal approach to single-observation bandwidth synthesis, as developed 
by A.E.B. Rogers [4], is based on coherent addition of the fringes from the 
recorded channels within an RF band. In contrast, the BLRO approach is based on 
the separate extraction of phase for each channel and the e^liclt manipulation 
of those phase values to obtain BUS delay. This single-channel-extraction approach 
can be carried out without phase calibration [6], which was a necessary capability 
in the early stages of BLRO development. 

The BLRO single-observation approach relies on two conditions. First, except 
for cycle ambiguities, the BSA delay and all chamnl-pair BUS delays must be 
approximately equivalent in expectation value. (Compare Eq. 74 and Eq. 81b.) 

Second system noise error in channel phase must be much less than <me cycle (in 
practice, 0.03 cycle). In this appendix, we will neglect charged particle 
effects and will also asstsne that dispersive effects before the tone-injection 
point can be calibrated. Although charged particle effects violate the first 
condition in the strictest sense, the condition is adequately met for spanned 
bandwidths of the order of 100 MHz or less at S-band. 

When these conditions ate met, single-observation bandwidth-synthesis pro- 
ceeds as follows. After extracting the phase for each channel, one computes a 
composite BWS delay for each channel pair from the difference phase between 
channels as described in Section XII. The confuted delays will possess inte^r- 
cycle aotbiguities as represented by the last term in Eq. 81b. These ambiguities 
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In delay are removed by iteratively proceeding from the unambiguous but imprecise 
BSA delay to the BHS delay of the most widely spaced channel*^alr» using at each 
step the result generated by the preceding channel pair. ' At each step, one com- 
putes the number of cycles that must be added to the BUS delay to place it 
within ±0.5 cycle/ of the (k-l)*-" delay, where Affe is the frequency separation 
for the kth channel pair. 

As an illustration of the bandwidth s}nithesis process. Table HI presents 
for our example observation (see Figs. 3, 5 and 6) the BSA delay and the three 
BUS delays obtained at five stages in the phase calibration and bandwidth synthe- 
sis process: (a) the uncorrected delays derived from the fringes (PHASOR output), 

(b) and (c) the delay corrections derived from tone phase (PCAL output) for each 
antenna, (d) the phase-calibrated but ambiguous delays, and (e) the final unambig- 
uous delays. For simplicity, only the BSA delay from channel 1 is used, although 
averaging all three S-band channels should improve the BSA delay determination. 

The tone BSA delay for each antenna was obtained from the phase-frequency slope 
of tone phase, which was equal to -11.221 psec for the example of DSS 43 in 
Fig. G2. An example calculation of tone BUS delay was given in the preceding 
appendix for channels 1 and 3 for DSS 43. Uhen the delay corrections are com- 
bined with the uncorrected delays, one obtains delays that are Ideally free of 
unwanted Instrumental effects but possess integer-cycle ambiguities. As dis- 
cussed above and as schematically Indicated in Fig. HI for a hypothetical obser- 
vation with low SNR, single-observation bandwidth synthesis can now be applied 
by iteratively removing ambiguities while progressing from the leftmost to the 
rightmost value for delay in Table HI. The final value for each of the four 
delay observables should be the residual delay given theoretically by the com- 
bination of delays in Eq. (77). Note that, while the final delay values for 
the various channel pairs in Table HI agree fairly well, they are not in agree- 
ment at the system noise level, as one might hope. This disagreement is probably 
a consequence of one or more of the following dispersive effects: (a) Improperly 

calibrated nonlinear trends in the phase-frequency response of the instrumenta- 
tion (see preceding appendix), (b) dispersive effects in the input tones and (c) 
dispersive effects such as multipathlng before the tone injection point. It is 
outside the scope of this report to estimate such errors. However, it should be 
pointed out that the differences between BUS delays are not necessarily an indi- 
cation of the error that will be propagated to geophysical/astrometric param- 
eters since large components of the differences would be constant and those 
components would therefore be absorbed by the clock parameter during multi- 
parameter estimation. 

A parameterized error analysis of the above BUS process leads to an asso- 
ciated method for computing a channel placement that helps in minimizing the 
nunri>er of inner channels and thereby Increases precision. In this particular 
method, it is assumed that a priori information is too poor to help with any 
step in the ambiguity resolution process (except tone phase slope). To reliably 
compute integer cycles for the (k+1)*^^ channel pair, the error ( 0 |j) in the delay 
for the k^*' channel pair and the frequency separation (Afjj+i) for the (k+l)^h 
channel pair must satisfy the relation 

M ; 0,5_acli (Bl) 

k+1 k 3 
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KAOeUY WnOElAY MADEUY 

NOAMMGUiTiS NO. 1 N0.2 

Fig. HI. Schematic example of ambiguity resolution In 
the BLRO system 


103 


Table HI. An example of phase calibration and single-observation bandwidth 
synthesis with the present BLOCK 0 syston 
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Relative to central lag of correlator in nanoseconds. Errors are system noise errors. Channel refer- 
ence frequencies are (f , , f., f.) - (2276, 2299.2, 2271) MHz. 



where the factor of 3 serve« to Insure a 3o probability of success. The error 
in delay for the pair is given by 


/ 2 

[2a + 4of) 
(2ir Af ) 


(H2) 


where Og is the system noise error in phase (In radians) for a single channel 
(but involving two stations) and (radians) is equal t.' the RSS of all other 
dispersive errors in phase for a single channel from one station. The term 
2o| accounts for two channels while 4a| accounts for two channels at two 
stations. For simplicity, we assume that dispersive errors have the same magni- 
tude for both stations and all channels and are uncorrelated between channels. 
The system noise error in phase (radians) for a single-sideband system is given 
by 


a 

s 


1 

vT SNRg 


(H3) 


wher£ SNRg is the SNR at peak fringe amplitude in the given channel. The factor 
of /2 fringe is the improvement in phase precision resulting from the use of all 
lags for the channel. 

Inserting the last two equations into Eq. (Hi), one obtains an expression 
for the maximum allowed ratio (R) of channel separations (frequency multiplier): 


k+1 

k 


^ R 


(H4) 


where 

j, = 1 y 

* ■ 3 T72 (H5) 

and where is still in radians. 

Given the single channel SNR and the other dispersive errors in phase, one 
can easily ccmtpute the frequency multiplier. Given the multiplier and the sepa- 
ration of the first channel pair, one can compute the successive channel separa- 
tions. However, the actual spread can increase somewhat more rapidly than the 
above multiplier implies. As illustrated in Fig. H2, pair 2 ^comprised of chan- 
nels 2 and 3) is given a separation of Af 2 ■ RAfj^ where Af| is the separation 
of pair 1 (channels 1 and 2). The separation of pair 3 (channels 1 and 3) is then 
given by Af 3 > (R-fl)Af]^. When a fourth channel is added, pair #4 (channels 3 
and 4) is given a separation of R times the largest previous separation or 
Af 4 ■ RAf 3 ■ R(R+l)Af 3 . The separation of pair 5 (channels 1 and 4) then becomes 
Afe ■ (R-fl)^Af,. Continuing this process, one finds that the outer pair separation 

^ X H ^2 

for a system with N^ channels would be given by (R't'l) ^ Af^. 
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(R-» 



Af^-RAf3>R(R + l)Af, 







CHANNEL NUMKR CHANNEL mEQU04CY ■ - » 

Fig. H2. An example of the geometric progression approach to 
BWS channel placement 
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This leaves us with the problem of specifying the maxi m um allovad 
separation of pair 1 In terms of the channel bandwidth. The delays obtained 
from bitstream allgmiient are based on the continuous channel bandpass and are 
therefore unambiguous. These delays can be used to Initiate the BW8 process 
by resolving the cycle ambiguity In the first channel pair, provided the 
separation of the first pair Is not too large. Without presenting all of the 
details, we will state here that a rough estimate for the multiplier from 
channel bandwidth to the first channel pair Is given by Afi “ R/1.8, where R Is 
computed by Eq. (H5). (The factor of 1.8 comes from (a) the fact that the 
continuous channel bandpass makes inefficient use of bits relative to the case 
«diere all bits are at the edges of the bandpass, and (b) the fact that tones 
placed within the channel bandpass for intrapassband phase calibration would 
not be placed at the bandpass edges but would be displaced inward by a safe 
amount. ) 

Combining the above results, one obtains an expression for the separation 
(Afn) of the most widely spaced pair in terms of the channel bandwidth (W) and 
the number (Me) of channels: 


W 


R 

1.8 


(IW-l)"c 


-2 


(H6) 


Relative to the first reference channel, the channel frequencies In units of 
the channel bandwidth would be 


n _R_ _R (R+1) _R_ (R+1)^ _R_ (R+1) ® (H7) 

1 . 8 ’ 1.8 ’ 1.8 ’ ” 1.8 

For convenience, the curves connecting Af)f(N(.), R and are plotted in Fig. H3. 
The right-hand vertical axis gives the total allowed dispersive phase error 
(numerator in Eq. H2) associated with each frequency multiplier R on the 
left-hand vertical axis. When using the graph to compute the required number 
of channels for a given spanned bandwidth and given frequency multiplier, one 
must adopt the next higher value of N,., if the point falls between two curves 
of constant N,,. 

The BWS approach outlined above is not optical for all applications. For 
smaller values of the frequency multiplier (R < 5), the basic assumptions behind 
Eq. H6 (e.g., single-channel detection) can be violated. Further, if good 
a priori information for delay is available, it may be possible to omit some or 
all of the inner channels required for ambiguity resolution. More generally, 
totally different techniques can be used in postcorrelation data reduction. For 
example, if detection and estimation were carried out through colwrent addition 
of BWS channels (Ref. 4), then one might choose to place the channels across the 
RF band in a pattern based on detection capability in addition to final delay 
precision. 
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Fig. H3. The relationship betweoi spanned bandwidth, BUS 
frequency multiplier and number of channels 
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APPESDIX I 


GEOMETRIC TIME DBUY 


'nils oppmdlx derives an .xpressloa for tte geoaetrlc delay in cetM of 
station locations and source direction. 

Suppose two antennas slnultanaously receive tlie slpial fron a very distant 
point aource. The dlffer«ace In the arrival tines at tha two stations Is 
referred to as the geomtric delay. Due to the great distance of the source, 
the signal can be represented as the plane wave In Eq. (7). All quantities will 
be refermced to geocentric coordinates In a nonrelativlstlc analysis. 

Let a particular wavefront rea^h station 1 at true tine t^. This segnent 
will reach station j at true time tj when the following phase equality Is satis* 
fled. 


k*x^(tp - oit^ ■ k*x^(tj) - «tj (II) 

Let tj • t^ -f Tg, where t. is the geometric delay. Rearranging Eq. II, we 
obtain 


o>T^ - k • tx^(tJ+Tg)- Xj(tp] (12) 

This equation defines Tg In terms of station location and source direction but 
It Is not computationally convenient. An accurate approximation can be used 
to obtain a more tractable f^rm. Since Tg is small (^.02 sec for earth*flxed 
antennas), expand xj about tj to obtain 


X (t!) - X (tJ) + X. T 

j J j 1 j S 


One can now use Eq. 


(12) to solve for t : 

O 



( 13 ) 


( 14 ) 


where f is the baseline vector given by 


B(tJ) - Xj(tJ) - x^(tp (15) 

As one irauld expect, the g«»etrlc delay Is sla^ly the Instmtttiainis path 
difference divided by the speed of light with a small correctlMi for the rnotlmi 
of station j during wavs transit. As noted in the text, this definition of the 
ge(»etric del^y uses station 1 as the reference statloa. That Is, If a s^Mot 


lOf 


of tho twye •trtkfo ototion 1 at tine ti« ehm that a^^mt will roach atatlon J 
at tine t^ * Tg(t^). where the tine tag for Tg is the arrival tlae at atatlon 1. 

A coonent concerning the neanlng of antenna position la In order. It 
is custMwry to make the "Intersection" of the two axes of rotation the point 
that defines the location of an antenna. Bven though the signal Is not received 
at that point, corrections can be applied to the observed delays to nake this 
point the effective point of reception. Two types of correction are require. 

The first is a correction applied only to those entennas for idilch the axes do 
not, in fact, intersect, as in a DSN equatorial nount. This correction d^ends 
on pointing direction, and it, in effect, shifts one of tiM axes of rotation in 
order to define a fictitious but earth-fixed point of intersection. The seco nd 
correction is a constat delay correction to account for the fact that the sig- 
nal is actually received at the feed rather than at the intersection of axes. 

This correction is equal to the difference betvraen the transit tistt that a wave 
would experience in its actual propagation to the injection point and its 
theoretical transit tine to the intersection of axes. This correction is noninally 
the same for all sources since the antenna always pcints toward the source and 
therefore always presents the same geometry in the propagation direction. For 
very accurate measurements, this correction nay not be constant but may dep«^ on 
pointing direction. 
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